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ABSTRACT 


An  automated  cloud  analysis  program  was  developed  and 
established  on  the  SPADS  computer  system  at  the  Naval  Envi- 
ronmental Prediction  Research  Facility  (NEPR?) .  The  pro- 
gram evaluates  GOES  visual  and  infrared  satellite  imagery 
simultaneously.  The  analysis  method  produces  information  on 
cloud  types,  cloud  amount,  precipitation  intensity,  and 
cloud  top  height  and  temperature  through  use  of  threshold 
tests  of  radiance,  texture,  and  temperature.  A  review  of 
current  work  on  the  evaluation  of  satellite  information  by 
computer    and   by   manual  analysis    is  included. 

A  maritime  region  460  X  463  nautical  miles  in  size  was 
selected  for  test  analysis.  The  satellite  imagery  was  manu- 
ally evaluated  and  compared  to  the  computsr  generated  out- 
put. Reasonably  good  patterns  of  cloud  types,  precipitation 
and  cloud  amount  were  produced  by  the  computer,  although 
further   testing   and   verification    is   needed. 
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I-       INTRODUCTION 

The  personnel  operating  Naval  ships,  aircraft  and  weapon 
systems  need  near  real-time  sub  synoptic  and  mesoscale  mari- 
time weather  information  to  operate  effectively  and  safely. 
Current  conventional  observational  data  (from  surface  ,  raw- 
insonde,  and  aircraft  observations)  do  not  adequately  cover 
the  maritime  region  cr  provide  for  mesoscale  analysis  of 
weather  phenomena.  Specifically  needed  are  detection  and 
monitoring   of      any   adverse     weather   conditions.  This   need 

pertains  to  phenomena  ranging  from  mesoscale  (10-100  km)  to 
synoptic  scale  (100-1000  km)  which  are  poorly  resolved  by 
conventional  observations.  Satellite  observations  can  solve 
this  problem  through  the  acguisition  of  global  high  density 
data  which  are  available  regardless  of  local  influences  of 
geography,    surface   conditions,    or   local   politics. 

Another  problem  is  the  need  to  provide  an  analyzed  prod- 
uct to  the  user  in  a  short  (1-3  hour)  time  frame.  Manual 
analysis  of  satellite  data  is  time  consuming  and  tedious, 
and  it  does  nor  provide  the  information  in  a  real-time 
fashion.  An  interactive  computer  system,  on  the  other  hand, 
can  collect,  analyze,  and  produce  the  required  products 
within  the   required   time    limit. 


14 


This  research  effort  concentrated  on  the  development:  of 
techniques  to  specifycloud  types,  cloud  amount,  cloud 
height,  and  precipitation  intensity  with  an  emphasis  on  the 
maritime  regions  from  satellite  images.  Information  on 
critical  weather  parameters  such  as  low  visibilities,  ceil- 
ings, precipitation  presence,  and  intensity  can  be  derived 
from  these  cloud  characteristics.  This  work  uses  Geosta- 
tionary Operational  Environmental  Satellite  (GOES)  Visual- 
Infrared  Spin  Scanned  Radiometer  (VISSR)  digital  satellite 
data  from  visual  and  infrared  channels  which  have  a  one-half 
hour  temporal  and  0. 5  to  4  nautical  miles  (n  mi.)  spatial 
resolution.  The  data  were  received  and  processed  at  the 
Naval  Environmental  Prediction  Research  Facility  (NEPRF)  in 
Monterey,  California  using  the  Satellite  Data  Processing  and 
Display  System  (SPADS).  This  interactive  computer  system 
was  developed  at  NEPEF  for  use  by  the  operational  Navy  to 
give  real-time  information  about  weather  phenomena.  One  sys- 
tem is  now  operational  at  the  Naval  Eastern  Oceanography 
Center    (NEOC)     at   Norfolk,    Virginia. 

The  objective  of  this  thesis  is  to  establish  and  evalu- 
ate a  computer  program  on  the  SPADS  to  analyze  GOES  images 
to  derive  cloud  types,  precipitation  intensity,  cloud 
height,    and   cloud    amounts.        To    accomplish    this,    an    in-depth 
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study  of  previous  research  and  analysis  techniques  was 
initiated,  and,  from  this  study,  useful  algorithms  and 
techniques  were  gleaned.  These  algorithms  and  techniques 
were  combined  to  produce  the  desired  computer  program.  Pre- 
liminary tests  of  the  computer  program  were  completed  on  the 
SPADS   system    at    NEPR?. 

Chapter  II  begins  with  a  review  of  previous  studies. 
This  includes  discussions  of  four  techniques  of  cloud  type 
classification;  (1)  cloud  index  (Harris  and  Barrett , 1975, 
1978),  (2)  spectral  analysis  (Lil jas, 198 1  a,  1981b),  (3) 
3-dimens ionai  nephanalysis  (Fye,  1978),  and  (&)  two-dimen- 
sional histogram  (Piatt , 1 981)  .  A  discussion  of  the  techni- 
ques  for  determining  cloud  height  and  clsui  amounts  used  by 
Reynolds  and  Vonder  Haar  (1977)  and  others,  and  of  proce- 
dures for  identifying  precipitation  areas  and  their  associ- 
ated intensities  is  also  covered.  Chapter  III  describes  the 
set  of  algorithms  selected  for  implementation  on  SPADS. 
Chapter  17  presents  an  outline  of  the  test  procedures  used 
and  reports  the  results  of  these  tests  (including  detected 
problem  areas).  Chapter  V  summarizes  the  research  efforts 
and  problems  and  provides  suggestions  for  future  changes  and 
additional    research. 
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II.       CLOUD  CLASSIFICATION    TECHNIQUES 

A.       INTRODUCTION 

Satellite  visual  (VIS)  and  infrared  (IR)  imagery  can 
provide  a  wide  variety  of  weather  information,  including 
estimates  of  cloud  type,  cloud  height,  cloud  amount,  and 
precipitation  intensity.  In  fact,  a  satellite  nephanalysis 
can  produce  a  spatial  summary  of  these  important  cloud  fea- 
tures. An  analyst  examines  the  visual  brightness  to  derive 
estimates  of  thickness  of  the  cloud  and  the  cloud  texture  to 
differentiate  between  cumulus  and  stratus  type  clouds.  3y 
using  the  IR  imagery,  cloud  top  heights  can  be  estimated; 
the  brighter  the  cloud  image,  the  colder;  therefore,  the 
higher  it  is.  With  combined  visual  and  infrared,  one  has  a 
greater        ability  to        pick        out      the      cloud     types      and, 

furthermore,    to   estimate   the    probability  of    precipitation. 

One  of  the  main  problems  encountered  in  the  1960'  s  and 
early  1970's  was  the  lack  of  standardization  of  rules  gov- 
erning the  nephanalysis  construction.  Each  country,  and 
sometimes  each  analyst,  has  used  subjective  methods  for  pro- 
ducing cloud  analyses.  Some  of  the  specific  problems  aris- 
ing   from    lack   of  standardization,    as   discussed   in    Harris   and 
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Barrett  (1975) ,  have  been:  no  minimum  size  is  established  as 
reference  for  inclusion  or  exclusion  from  analysis;  cirrus 
clouds  are  rarely  identified  in  middle  latitudes; 
significant  cloud  areas  are  subjectively  delineated  at  the 
discretion  of  the  analyst;  only  four  cloud  categories  are 
assigned  (which  consisted  of  unegual  percants  of  areas)  ;  and 
no  standard  nephana  lysis  construction  rules  have  been 
published. 

With  these  problems  in  mind,  Harris  and  Barrett  (1975, 
1978),  Barrett  and  Martin  (198  1),  Piatt  (1981),  Fye  and 
Logan  (1977)  and  Lil  jas  (1981a)  ,  to  name  a  few,  have  been 
working  toward  the  goal  of  producing  a  recommended  standard 
for  objective  nephana lysis  of  satellite  VIS  and  IR  imagery 
through  manual  and  computer  analysis.  The  following  sections 
review  these  representative  research  efforts  in  this 
endeavor.  The  discussions  are  divided  into  three  areas  of 
concern  to  this  thesis;  cloud  typing,  cloud  height  and 
amounts,    and   precipitation. 

B.       CLOUD    TYPING 

Basically  four  distinct  approaches  to  effective  cloud 
typing  have  been  proposed  in  the  literature.  They  are:  (1) 
cloud   index    (Harris    and    Barrett,         1975,       1978;       Barrett    and 
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Harris,  1977)  ,     (2)    the   O.S.    Air  Fores    3DNEPH    program    (Fye   and 
Logan, 1977;        Fye, 1978),         (3)         spectral      cross-correlation 
(Liljas, 1981a,       1981b),      and      (4)      two-dimensional   histogram 
(Piatt,  1  98  1)  .      These    methods  are    reviewed   below. 
1 .      Cloud    Index 

Harris   and      Barrett    (197  5)      recognized  the     need   for 
standardization    and   fcr    increased   information   content.       As    a 
result,      they   developed      a  cloud    indexing   system      which    pro- 
vided specific      guidance   and      rules   for     analyzing   satellite 
visual  imagery      manually.        For      clarity   and      understanding, 
they    proposed      the   analysis      be    composed      of   three      separate 
layers  of      analyses;      two    descriptive   and      one    interpretive. 
Fig.      1    is    the   flow   chart     of   the  stages   in    the   construction 
of   this   nephanalysis.      The   authors  developed    a    key    (Fig.      2) 
which   includes   details  on   the   following   cloud   features;       (1) 
percent   of    cloud,  (2)     cloud  type     (   e.g.      cumulus),     (3)       cloud 
structure    (including      size,      shape,      and   pattern),        and    (4) 
heighx  of      cloud        tops.         The      key  also      gives    interpretive 
weather   phenomena. 

The  first  analysis  layer  contains  information  on 
geography  and  cloud  cover.  The  cloud  cover  is  divided  into 
five  egual  percent  categories  with  a  minimum  size  reported 
of    a    2   1/2    degree   sguare. 
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The  second  descriptive  lays r  includes  cloud  type 
and  cloud  structure.  They  are  classified  in  the  following 
order:  major  cloud  systems,  definite  boundaries  and  finally 
indefinite  boundaries.  The  analyst  uses  cloud  brightness 
and    texture  to   discern  six  types    of   clouds    (see   Fig.    2) . 

The  final  analysis  layer  produced  is  the  interpreta- 
tion of  the  cloud  features  using  the  information  from  the 
first   two    layers.  The   types   of  cloud      features  considered 

important  are  displayed  in  Fig.  2.  All  three  analysis  layers 
combined  were  designed  to  be  completed  within  one  hour  by  a 
skilled  analyst. 

Barrett  and  Harris  (1977)  extended  this  procedure  to 
include  infrared  imagery  for  the  purpose  of  supplementing 
the  visual  analysis.  This  added  capabilties  for  nighttime 
and  high  latitude  analyses.  They  mention  three  possibilities 
for  analyzing  IR  imagery;  (1)  use  longstanding  internation- 
ally accepted  codes  and  symbols,  (2)  follow  the  visual  pro- 
cedures previously  outlined,  or  (3)  use  a  new  procedure 
which  recognizes  and  represents  the  special  properties  of 
the  13  images.  The  first  two  methods  would  be  directly  com- 
parable to  the  visual  analysis.  The  third  approach  would 
offer     additional      information   to      allow      identification      of 
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Satellite  image 


Outline  and  denote  areas  of 
uniform  %    cloud  cover  in  the  order  :- 

a)   Major  cloud  systems 

b)  Cloud  areas  with  detinue  boundanes 

c)  Cloud  areas  with  indefinite  boundaries 
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Satellite  nephanalysis 


Figure   1. 


Manual    Nephanalysis 
Barretx,    1 975) . 


lew   Chart    (from    Harris    and 
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Figure  2.      VIS   Key  to   be    Dsed   in    Manual   Nephanalysis    (from 
Harris   and   Barrett,    1975). 


cloud  height  through  cloud  brightness  evaluation  and  a  more 
specific   identification   of  cloud    texture. 

The   IB      analysis    technique   identifies      the   following 
cloud   features    (see  Fig.      3):  (1)       percent   of   cloud   cover 

(same  as  VIS  technique),  (2)  brightness  of  cloud,  (3)  cloud 
texture,  and  (4)  size  of  clouds  (same  as  VIS)  ,  plus  inter- 
pretive categories  with  eighteen  of  then  the  same  as  VIS 
(heights  and  three  cloud  descriptions  ire  not  included). 
Barrett  and  Harris  (1977)  felt  that  any  of  the  three  methods 
could  produce  an  analysis  in  a  single  hour  if  done  by  a 
trained   analyst. 
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Figure  3.      IE    Key  to    be   Used   in    Manual   Nephanalysis    (from 
Barrett   and  Harris,    1977). 


After  establishing  a  manual  method  for  both  VIS  and 
IS  satellite  nephanal  ysis,  Harris  and  Barrett  (1978)  devel- 
oped an  automated  objective  nephanalysis  through  a  cloud 
rscognition   scheme.  They   used    an   approach      which   axamined 

guantitative  cloud  brightness  and  texture  within  a  small 
area  of  the  whole  picture.  The  flow  chart  in  Fig.  4  out- 
lines their   approach.      Cloud    brightness   was    assessed   in    each 
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subarray  by  summing  the  brightness  counts  of  only  cloud  pix- 
els and  then  dividing  by  the  number  of  cloud  pixels  (Eq.  1 
and  2,  Table- I).  Cloud  amount  was  calculated  in  each  subar- 
ray by  summing  the  number  of  cloud  pixels  and  then  dividing 
by  the  total  number  of  pixels  per  subarray.  This  number  was 
then  multiplied  by  100  to  obtain  a  percent  of  cloud  cover 
(Eg.    3  in    Table   I)  . 

Their  assessments  of  cloud  texture  used  statistical, 
instead  of  structural,  measures  of  local  variation  in  image 
density.  One  measure  they  employed  was  the  standard  devia- 
tion of  each  subarray  (5x5  pixel  array)  of  density  values, 
as    calculated    with   Eq.   4     (Table    I). 

Their  other  criterion  was  vector  dispersion,  which 
required  that  the  density  values  be  trsated  as  a  set  of 
adjacent  triangular  planes.  The  dispersion  in  three  dimen- 
sions of  the  normal  tc  these  planes  gave  another  measure  of 
texture.  Fig.  5  shews  a  graphical  representation  of  -his 
method.  Eq.  5,  6,  and  7  in  Table  I  were  used  to  calculate 
these   estimates    of   texture. 

Once  all  these  values  were  calculated,  boundaries 
between  cloud  types  had  to  be  established.  Three  categories 
were    selected;       (1)       sheet   or    layered      clDuds      (stratiform), 
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(2)  cellular /tower  clouds  (cumulus)  and  (3)  broken/mottled 
cloud  (s tr at o cumulus  or  mixed)  (Harris  and  Barrett, 1978) .  A 
three-dimensional  decision  spaca  was  sat  up  using  two  tex- 
ture measurements  and  one  brightness  vaiue.  Those  areas  that 
had  a  brightness  above  the  threshold  were  included  in  fur- 
ther analyses;  those  below  were  not.  Fig.  6  shows  the  deci- 
sion space  based  on  the  two  texture  measurements.  The 
discriminant  lines  were  developed  through  analysis  of  a 
training  set  of  Defense  Military  Satellite  Program  (DMSP) 
sample  brightness  and  texture  measurements  of  each  cloud 
category . 

In  assessing  the  accuracy  of  this  method,  Harris  and 
Barrett  (1978)  found  several  problem  areas.  For  example, 
towering  cumulus  in  post  frontal  areas  (especially)  were  not 
identified  because  of  the  field  of  view  (FOV)  size  and  snow 
areas  were  included  as  clouds.  But  the  overall  accuracy  was 
still  good  with  a  greater  than  12%  correct  classification. 
Harris  and  Barrett  (1978)  felt  that  several  procedures  would 
improve  the  output  of  the  model.  Among  these  were  increasing 
the  resolution  of  the  data  by  decreasing  the  size  of  scan- 
ning spot,  decreasing  the  subarray  size,  normalizing  the 
image     brightness    for      changes   in      sun      viewing   angle,         and 
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increasing  the     texture    parameters  to     permit   identification 
of    more   cloud  classes. 
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Figure  4.      Flow  Chart    for   Automated  Objective    Nephanalysis 
(from   Barrett    and   Hams,     1378)  . 

2-      Three-Dimensional    Nephanalysis    (3DNEPH) 

The   United      States    (CJ.S.  )  Air    Force      developed    an 

automated  three-dimensional  cloud  analysis  called  the 
"3DNEPH  model".  It  uses  satellite  imagery  from  the  DMS?  and 
TIROS-N   satellite   series    together     with   conventional   in   situ 
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a) 


V 


Figure  5 . 


Graph   of   Vector   Dispersion   Techniaue.       (a) 
represents  a    smooth* s urface    (stratiform   clouds) 
ana    (b)    represents   a    rough    surface    (cumulus 
clouds).    Normal   vectors   ror   each    case  are 
included   to  the  left    of    the   diagrams     (from    Harris 
and    Barrett,    1978)  . 


data  from  aircraft,  surface,  and  upper  air  observations.  A 
cloud  analysis  is  produced  in  a  reai-tima  fashion  regardless 
of  ice,  snow,  desert,  or  other  geographical  anomalies  (Fye 
and  Logan, 1977) .  The  satellitas  provida  high  density  data 
over  tha  entire  world,  thus  permitting  a  global  analysis  of 
clouds.  The  horizontal  resolution  of  the  output  is  25  n  mi. 
The  data  array  has  a  pixel  resolution  of  3  n  mi. 
The  sataliite  data  are  first  formatted,  rectified,  mapped, 
and  stored  on  a  computer  storage  device  before  running  the 
3DNEPH  program.        The   rectification   corrects    for   the   effects 
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TABLE    I 

Aloqrithms   for   the   Objective  Nephanalysis    (from   Harris   and 

Barrett,     1978) 


(1)  N  N 

„      1-1    j-1 

B   ■  •> 


(2)  J!  N 

i-i  j=i 


B  -  cloud  brightness 

N  -  number  of  rows/columns 

X..  -  density  at  each  point 

6..  -  step  function  (  -  1  X,  _>  t 


ij 


;  -  o  %..  <  t 


ij 


t   -  density  threshold  defining  the 
boundary  between  cloud/no  cioud 


(3)  A  =■  100  (M/N") 


A   -  average  cloud  amount 


(4) 


N    N 


a  -  [T^  ]|P     CXtj  -  x)2/N2  -  l]1/2 


i-1    j«l 


a   -  standard  deviation  of  each  sub- 
array  density  values 

X   ■  mean  of  the  subarray  density  values 


«   t-iat 


(£.,  m  ,  n  )  ■  direction  cosines  of  the  ith  plane 


*  "  R  Eni 


(£,  m,  n)    -  max  likelihood  of  true  polar 
vector  (X,  y,  y) 


n  -  -  Lni 


(6)  R2  -  (IJ^)2  +  (Im^2  +  (En^2 


(7)   k 


N-l 

N-R 


where: 


N   -  the  sample  size 

k   =  estimate  of  the  texture  of  the  surface 

k-1  =  smooch  surface  (stratiform) 

k~0  =  uneven  surface  (cellular  cumulus) 
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Figure  6.      Decision   Space    Using    Two   Texture    Measurements, 
Standard  Deviation   and   Vector   Dispersion.    Line 
('I   separates    stratiform   from   mixed   statccumulus 
ana   line    (2)     separates    mixed    from   cumulus    clouds 
(from   Harris   and  Barrett,    1973)  . 
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of  limb  darkening  before  storage.  Limb  darkening  is  the 
effect  of  a  decrease  in  brightness  as  the  edges  or  "limbs" 
of  the  earth  scan  are  approached.  The  data  are  not  normal- 
ized with  respect  to  light  (VIS)  or  longwave  radiation  (IR) 
(Fye  and  Logan, 1977).  Each  visual  and  infrared  pixel  is 
represented  by  a  single  grayshaie  value  in  the  range  from  1 
to  63.  In  the  visual,  1  is  dark  and  63  is  white.  Infrared 
pixels  represent  temperature  in  degrees  Kelvin  (K)  from  210 
K  for  a  value  of  1  to  310  K  fac  63;  eaoh  grayshade  repre- 
sents a  1.6  degree  change  (Fye,  1978)  .  The  preliminary  pro- 
cessor calculates  the  average  grayshade  value  per  25  n  mi. 
square  using  Eq.  1  (Table  II)  and  then  uses  this  to  get  the 
variability  within  the  square  using  Eq.  2  (Table  II) 
(Fye, 1978).  Both  these  calculations  are  performed  sepa- 
rately for  13  and  VIS  images.  The  VIS  and  IR  grayshade  val- 
ues and  variabilities  are  used  to  identify  the  various  cloud 
types.  Fig.  7  and  8  represent  the  cloui  typing  algorithms 
used.  They  were  constructed  empirically  by  comparing  many 
grayshades  and  variabilities  to  corresponding  satellite 
imagery  and  surface  data  (Fye, 1978).  The  operational  algor- 
ithms were  derived  from  the  figure  and  then  modified  to 
obtain  optimum  results.    One   obvious  problem  with  electing 
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to  use  the  model  (given  the  available  documentation)  is  the 
lack  of  graduated  values  on  the  axes  of  these  two  figures. 
Only  relative  values  of  warm/cold  or  light/dark  are  included 
in  the  report.  Also  lacking  in  the  Fye(1978)  report  are  any 
details   of   the   current  algorithms   used   for    cloud   typing. 


TABLE    II 

3DNEPH   Statistical    Eauations    for   the   Average   Grayshade   and 

Variance    (from    Fye,    1978) 


N 
(1)      G   =   1/N    £      G. 
i=l 


N 


,   G>0,    N>0 


(2)      V  -  1/N  Y*    |G     -  G|    ,  G>0,   N>0 
,      i=l 


G  =  average  gray  shade 

G.    =    individual   3   nmi  area 
grid   gray   shade 

N  =    total   number   of  G.    per 
25  nmi  grid   space 


V  =  variability  within   the 
25   nmi  grid    space 


-•      Spectral  Cross-correlati  on 

The  third  classification  method  is  the  spectral 
cross-correlation  method  developed  by  Erik  Liljas  (1981a) 
for  use  on  an  interactive  computer.  Liljas  used  TIROS-N  and 
National  Oceanic  and  Atmospheric  Administraion  (NOAA)  6 
satellite  data  from  three  spectral  channels;  channel  1 
(VIS)  ,  channel  2  (near-IB)  and  channel  'i  (IR)  .  During  the 
summer  time,    different  cloud   types,      land    and    water   surfaces 
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Figure  7.      The   Decision    Space   of    the    3DNEPH   Program.      IR   and 
and   VIS   grayshades   are   used   to    determine  the 
cloud   type    (from   Fys,     1978). 
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Figure   8.      An    Additional   Decision   Space   of    the    3DNEPH 
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cloud   typing    (from    ?ye,    1978). 
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were  analyzed  to  detect  their  unique  spectral  signatures. 
The  results  were  used  to  build  a  three  dimensional 
parallelepipedon  classification.    Fig.    9. 

Two  satellite  types  were  used  far  the  test  of  the 
model;  polar  orbiting  (TIROS-N  and  NOAA-6}  and  geostationary 
(Meteosat)  .  The  satellite  data  used  for  the  test  of  the 
model  were  geometrically  corrected,  but  they  did  not  need  to 
be  normalized  because  the  mid-day  pass  of  the  satellite  was 
used        (approximate        sun        elevation  of        45        degrees). 

Normalization  corrects  for  the  radiance  variation  caused  by 
the  sun  location  (elevation).  Also  the  amount  of  data  was 
reduced  from  1#02U  digital  levels  to  256  levels,  which  was 
considered   adequa.-e      resolution    for      the    study.  Data    from 

channel  2  were  primarily  used  to  separate  land  and  waxer 
areas.  Channels  1  and  4  were  the  principal  cloud  *ype 
identifiers. 

Comparison  of  the  spectral  cross-correlation  method 
with  synoptic  observations  showed  that,  the  computer  method 
gave  correct  indication  of  cloud  types,  and  also  an  improved 
cloud  division.  As  shown  in  Fig.  10  (see  Table  III  for  sym- 
bol definitions  of  cloud  types)  ,  several  cloud  types  are  in 
the      same      boxes;         Liljas   recommends     the      use      of      texture 
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observation  to  discern  between  these  cloud  types.  This  was 
done  manually  and  entered  interactively.  Liljas  (1981b)  also 
reports  that  good  information  on  cloud  structure  of  weather 
systems  can  be  retrieved.  He  recommends  the  use  of  multiple 
parallelepipedon  classifications,  each  classification  asso- 
ciated with  a  different  sun  elevation  instead  of 
normalization   of   the    reflected   imagery. 
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'iaure  9.      Parallelpicedon  Classification    of. Cloud  Types   by 
IE   and   VIS "Digital   Counts    (rrom    Liljas,    1981a). 


34 


100 

VISUAL 


Figure   10. 


^^f^TTT^11^   Decision    Space    for   Typing   Clouds 
(Table   II  I. defines  t he    symbols    for    clouds)     and 
Precipitation    Intensity.       it   is    represented    by 
une   area    with    darK   diagonal   lines    (see   Tabl<=   x 

S«liie(fgS"!H5lidlfifiP*ion  of  the  i»t««ity 
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TABLE    III 

Cloud  Classification  to   be  Us=i    with   Fig.     10 

1981a) 


(from   Liljas, 


Main  and  Cloud  Types: 

1.  Cumulonimbus 

2.  Nimbostratus 

3.  Cirrostratus 


4.   Cumulus  conges tus 


5.   Stratocumulus 


6.  Haze/Stratus 


7.   Land 


8.   Water 


Storm  cloud  with  high  top 

Squall  cloud  with  scattered  showers 

Large  vertical  thickness 
Rather  low  topside 

Dense  cirrostratus 

C  irrus 

Thin  cirrus  over  water 

Dense  altostratus 

Thin  altostratus  over  water 

Thin  altostratus 

Dense  altocumulus 
Large  piled  up  cumulus 
Rather  small  and  flat  cumulus 
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U.      Two-dimensional    Histocrra  m 

Another  spectral  method  was  developed  by  Piatt 
(1981)  using  a  two-dimensional  histogram  representation  of 
cloud  types.  He  selected  satellite  data  with  typical  cloud 
fields  to  develop  a  histogram  decision  space  representation 
of  cloud  types.  Piatt  concentrated  on  evaluating  cloud  sys- 
tems instead  of  isolated  clouds  and  therefore  his  results 
are   biased   toward   synoptic  features. 

The  measured  visual  radiances  were  compared  with  the 
direct  radiance  from  the  solar  disk  on  the  satellite  to  pro- 
duce a  calculated  isotropic  albedo  per  2.5  km.  These  aibe- 
does  are  expressed  as  a  fraction  from  0  to  1.  The  infrared 
radiances  are  converted  tc  black  body  brightness  temperatures 
ranging  from  193.3  to  301.5  K.  A  histogram  was  then 
constructed  by  sorting  these  data  into  983  different  bins  of 
temperature  and  albedo.  Each  area  evaluated  was  approxi- 
mately 500  X  500  km,  with  roughly  40,000  values  per  histo- 
gram. The  total  bin  counts  were  plotted  on  a 
two-dimensional  grid  and  then  contoured  by  same  bin  count 
(see  Fig.  11).  This  cloud  histogram  representation  could 
then  be  used  to  classify  cloud  types  by  matching  plots  of 
the      satellite    data.  After    classifying      various   types      of 
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clouds,    Fig-       12   was    compiled   to    give   an   idealized  depiction 
of   the  location   of  clcud   types  by  albedo   and   temperature. 

Piatt  (1981)  discussed  several  factors  that  tend  to 
reduce  the  usefulness  of  Fig.  12.  One  cloud  problem  arises 
when  there  are  either  breaks  in  a  cloud  deck,  or  partially 
filled  field  of  view  (FOV)  in  the  satellite  view  area. 
Additionally,  clouds  of  the  same  optical  depth,  but  exhibit- 
ing breaks  will  have  a  different  albedo  measurement  than  an 
unbroken  cloud  of  the  same  optical  depth  (stratocumuius  in 
Fig.  14  demonstrates  this) .  Another  factor  that  changes  the 
albedo  and  temperature  values  is  variation  in  cloud  top 
height.  All  cf  these  factors  create  ambiguity  in  determining 
the  cloud  type,  amount,  optical  depth,  and  cloud  top 
temperature.  As  a  result,  the  author  recommends  that  addi- 
tional work  should  be  done  with  these  results  before  using 
the    two- dimensio  ral   histogram    approach    to    identify   clouds. 
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Figure   11.      Two-dimensional  Histogram    Plot    of    Albedo   Versus 
Temperature     (representing   cumulus    clouds)  (from 
Piatt,    1981)  . 
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Ideal  Histograms  of  Cloud  Types  (from  Piatt, 


1981). 
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C.       CLOUD    HEIGHT    AND    AMOUNT 

Thrse  methods  of  deriving  cloud  height  and  amount  have 
been  proposed  in  current  research.  Two  methods  are  an  exten- 
sion of  those  reviewed  in  the  cloud  typing  section;  cloud 
index  {Harris  and  Barrett,  1975,  1978;  Barrett  and  Harris, 
1977)  and  the  U.S.  Air  Force  3 DNEPH  (Fya  and  Logan,  1977; 
Fye,  1978).  In  addition,  a  new  bispectral  technigue 
(Reynolds   and   Vonder    Haar,    1977)     is  discussed. 

1.      Cloud   Index 

In  two  separate  studies,  Barrett  and  Harris  (1977; 
Harris  and  Barrett, 1 975)  developed  a  recommended  standard 
cloud  height  and  amount  determination  procedure  for  manual 
analysis  of  satellite  visual  and  infrared  images.  The  vis- 
ual method  combined  brightness,  cloud  type,  and  synoptic 
situation  to  determine  gualitative  estimates  of  cloud  height 
(low,  middle,  or  high  level).  Later,  they  developed  a  sys- 
tem for  identifying  cloud  levels  through  the  degree  of  cloud 
brightness  in  the  IB  image.  IS  cloud  information  specifi- 
cally lends  itself  to  this  task  since  the  radiation 
temperature  of  the  cloud  is  a  function  of  the  cloud-top 
altitude    (Barrett   and   Harris,    1977). 
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Some  problems  arose  with  cirrus  and  cloud  elements 
smaller  than  the  field  of  view.  In  this  instance,  satellite 
measured  temperatures  can  appear  warmer  than  actual  cloud 
temperature,  because  the  satellite  is  also  receiving  the 
radiation  temperature  of  the  surface  or  lower  clouds-  This 
problem  can  be  minimized  by  using  IB  and  VIS  together  to 
identify  cirrus  clouds  and  by  using  higher  resolution  data  , 
9.3.,  data  from  the  Advanced  Very  High  Resolution  Radiometer 
(AVHER)    on    NCAA    satellites. 

Fig.  3  is  the  key  for  the  IE  nephanalysis  and 
includes  the  cloud  brightness  steps  and  Fig.  13  is  an  exam- 
ple   of  how      these  are   applied.      This    analysis      can   be    inter- 

* 

preted  in  terms  of  gualitative  heights,  low  (dim)  _  to  high 
(very  bright)  (Barrett  and  3  arris, 1977)  .  cloud  amount 
determination  categories  are  also  delineated  in  Fig.  2  for 
Vis  and  Fig.  3  for  IR.  It  is  done  in  categories  of  size  of 
cloud  element,  e.g.,  100-150  km.  Harris  and  Barrett's 
(1978)  objective  automated  method  does  not  specifically 
address  cloud  height,  but  does  calculate  an  average  cloud 
amount  for  each  subarray  size  (5x5  pixels)  using  Eq.  3 
(Table  I)  . 
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2.      Three-Dimensicnal   Nephanalysis    (3DNEPH) 

The  U.S.      Air   Force      3DNEPH    program's    main  objective 

* 

is  to  produce  information  on  cloud  height  and  cloud  amount. 
It  uses  a  standard  vertical  grid  of  15  layers,  which  range 
from  the  surface  to  100  millibars  (mb|  .  The  IR  image  is 
used  to  determine  height  and  ela ment  cloud  amounts,  and  the 
visual  and  infrared  are  used  together  to  determine  total 
cloud  amount.  This  information  is  than  combined  with  sur- 
face, aircraft  and  rawinsonde  data  to  determine  the  final 
cloud   height    and   cloud  amount   analyses. 

The  preliminary  processor,  as  described  in  the  cloud 
type  section,  prepares  the  satellite  data  for  height  deter- 
mination. In  addition,  it  sorts,  screens,  and  combines  the 
conventional  data  received  from  the  Air  Force's  global 
weather  data  base  to  produce  a  vertical  and  horizontal 
3DNEPH  g ridded   layered  format. 

Next,  the  satellite  data  goes  to  the  visual  and 
infrared  processors.  The  VIS  processor  steps  are  outlined  in 
Fig.  14.  Each  pixel  grayshade  is  directly  compared  to  a 
grayshade  (background/surface  brightness)  specified  for  each 
25  n  mi.  square.  This  value  represents  a  clear  sky  earth 
radiance    for   each   square.      The   difference    between   each    pixel 


43 


grayshade  and  the  background  brightness  is  used  in  Fig.  15 
to  decide  if  there  is  a  cloud  or  not.  The  curve  in  Fig.  15 
represents  the  brightness  threshold  for  clouds;  any  vari- 
ability greater  than  this  is  a  cloud  area.  This  curve 
incorporates  the  background  brightness  and  the  varying 
brightness  associated  with  different  surfaces  (  e.g., 
desert).  The  total  number  of  cloud  decisions  divided  by  the 
total  number  of  available  pixels  represents  the  amount  of 
cloud  cover   in   the  square. 

The  infrared  processor  incorporat2s  a  slightly  dif- 
ferent analysis  method.  The  major  steps  of  the  processor 
are  included  in  Fig.  16  as  a  flow  chart.  The  first  step  is 
to  construct  a  frequency  distribution  of  the  6U  raw  gray- 
shades  for  each  25  n  mi.  area  to  detect  significant  modes 
(maxima)  of  6  or  mere  grayshades  or  12  or  more  adjacent 
shades.  If  more  then  2  modes  are  identified,  the  modes  are 
combined  on  the  basis  of  relative  amplitude  (strongest  modes 
are  retained)  and  the  differences  between  modes.  If  there 
is  no  apparent  mode,  then  the  processor  defaults  to  a  cloud 
threshold  temperature  determined  from  appropriate  surface 
temperatures.  3efore  any  cloud  decision  is  made,  a  correc- 
tion  is      applied    to      the    IR      temperatures    to      compensate    for 
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atmospheric  absorption  and  other  effects  such  as  zenith  and 
locaj.  angle  which  create  biases  in  the  IR  data.  Fig.  17  is 
an  example  of  a  curve  used  to  derive  bias  correction  for  a 
particular  satellite.  The  actual  bias  correction  is  the  dif- 
ference between  the  best  fit  curve  and  the  diagonal  line  of 
perfect  correlation.  Next  a  representative  temperature  is 
selected  for  each  mode  (usually  the  coiiest  in  the  mode)  and 
these  data  are  compared  to  a  gridded,  conventionally  derived 
surface  temperature  data  base  to  make  a  cloud/no  cloud  deci- 
sion. Fig.  18  graphically  represents  this  decision  process. 
The  final  cloud  decision  is  made  by  using  Fig.  19.  This 
compensates  for  irregularities  in  the  surface  temperature 
data  base  which  limit  the  reliability  of  the  cloud/no  cloud 
decision.  The  cloud  amount  is  determined  by  the  number  of 
pixels  in  the  cloudy  mode  (Fye  and  Logan, 1977) .  Addition- 
ally, the  13  processor  determines  the  cloud  height  by  com- 
paring the  coldest  temperature  in  the  cloudy  mode  to  a 
real-time   upper   air   data    base. 

Other  processors  include  the  conventional  data  pro- 
cessor, the  data  integration  processor,  and  the  manual  data 
processor.  The  conventional  data  processor  takes  the  sur- 
face,       aircraft,      and     upper      air      reports    and      produces      a 
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Figure   16. 


3DNEPH  Flew  Chart  of  the  IR  Satellite  Dara 
Processor   Program    (from   Fyef    1978). 
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horizontally  gridded,  eight  mesh  data  base  of  cloud  informa- 
tion for  15  vertical  layers.  These  layers  of  data  and  the 
output  of  VIS  and  IB  processors  are  merged  by  the  data  inte- 
gration processor.  The  processor  also  performs  meteorologi- 
cal consistency  checks.  The  cloud  amounts,  heights,  and 
types  from  VIS  and  IE  data  ar=  conver-sd  to  specific  cloud 
layer  amounts  and  heights.  Conflicts  ars  also  resolved  by 
this  processor  by  checking  age  of  data,  greatest  cloud 
cover,      layer      amount,      ana  cload      zop   height.  Next,      the 
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satellite  and  conventional  data  are  compared  and  integrated 
with  a  continuity  field  to  complete  the  cloud  analysis 
(Fye,1976).  The  continuity  field  is  the  previous  or  latest 
3DNEFH  analysis   which   fills  in  missing  data   and   parameters. 

The  analyst  next  enters  manually  encoded  input  to 
alter  the    3DNZPH      data  base  as   necessary.  The   manual    data 

processor  takes  this  information  and  incorporates  it  into 
the  overall  data  base.  This  provides  an  avenue  for  quality 
control,  correction  of  errors  and  improvement  of  the  overall 
analysis. 

3 .      Bispectral 

k  bispectral  method  for  determining  cloud  height  and 
amount  was  developed  by  Reynolds  and  Vender  Haar  (1977). 
Visual  and  infrared  data  were  aach  analyzed  quantitatively 
over  an  array  of  scan  spots.  This  analysis  yielded  cloud 
amount  and  cloud  top  temperature,  from  which  cloud  height 
was  derived  by  comparison  to  a  nearby  vertical  temperature 
sounding. 

The  method  was  based  on  manipulation  and  solution  of 
two  equations  (Eg.  1  and  2,  Table  IV) .  These  equations  are 
basically  comprised  of  adjustments  for  the  fraction  of 
cloud-covered  and  clear-sky  araas  in  a  single  field  of  view 
and    for      backbody   emissivity      of    the      clouds.         One      was   the 
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satellite   visual     radiance  equation    (Equation    1,        Table   IV) 
and    the   other      was  the  satellite    infrared      radiance   equation 
(Eq.    2,    Table   IV).    The  equations    related  the   amount   of   cloud 
(Acid   in   tenths)      to    either  cloud  albedo    (pcld   X   Hs)       in   the 
visual  or    cloud   radiances    (Ncld)  in   the    infrared  channels 

(note  that  Acid  plus  the  amount  of  clear  area,  Aclr  in 
tenths,  equaled  1.0).  The  IR  radiance  equation  also 
accounted  for  the  effect  of  emissivity  (e)  of  clouds  on 
cloud  and  surface  (Nclr)  radiances.  In  both  equations,  Aclr 
was  replaced  by  1.0-  Acid.  The  rest  of  the  terms  in  the 
equations  except  two  {Acid  and  Held)  were  available.  Averaqe 
values  of  the  cloud  and  surface  albedoes  were  known  and  more 
exact  values  can  be  obtained  in  practice  "from  the  actual 
array  of  satellite  measurements"  (Reynolds  and  Vonder  Haar, 
1977).  Measured  surface  temperatures  were  converted  to  IR 
radiances  and  used  for  the  surface  radiance.  The  emissivity 
was  initially  assumed  to  be  0.9  for  all  clouds.  This  left 
only  the  two  unknowns  to  be  determined.  2q.  1  (Table  IV)  was 
then  solved  for  Acid  and  Eq.  2  (Table  IV)  for  Ncld  radiance 
which  resulted  in  Eq.  3  and  4  (Table  IV).  The  Acid  value  was 
derived  first  from  the  visual  satellite  data.  Next,  Acid 
was    used    with   the   simultaneously      received    IR    data   to    obtain 
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the  cloud  radiance  value.  The  resulting  radiance  was  then 
converted  to  cloud  top  temperature  via  Plancks*  equation. 
The  final  step  was  to  use  a  local  upper  air  sounding  and  the 
cloud  temperature  to    obtain  the    cloud   top   height. 

Reynolds  and  Vonder  Haar  (1977)  performed  an  uncer- 
tainty analysis  on  this  technique  to  determine  on  what  the 
uncertainty  (A)  of  the  cloud  amount  (Acid)  and  cloud 
temperature  (zlNcld)  measurements  depend  and  how  large  an 
uncertainty  (A)  is  allowable.  It  was  shown  that  A  Acid 
depends  on  the  m  agnitude  and  precision  of  Hs,  Ms,  pcld,  and 
pclr,  where  4Hs  and  A  pclr  will  be  less  than  10%  of  the 
error.  Ncid  depends  on  the  magnitude  and  uncertainty  of  Ml, 
Acid,    Nclr   and    €. 

The  two  quantities  A  Acid  and  4Ncid  were  evaluated 
using  typical  mid-latitude  values  obtained  from  satellites. 
Table  V  gives  the  absolute  ^Acld  error  far  given  A  Hs,  Ms, 
and  pcld.  These  results  show  that  the  "uncertainties  must 
stay  below  the  10%  level  if  reasonable  accuracy  is  to  be 
maintained"    (Reynolds  and    Vonder    Haar,      1977).  This    accu- 

racy is  possible  through  use  of  ground  and  cloud  truth  meas- 
urements along  with  improvements  in  instrument  calibration. 
At    present      there   is    approximately      a    3    degree      Celcius    (°G) 
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TABL2    IV 

Algorithms   for   Cloud    Top   Temperature    and   3loud   Amount    (afte: 
Reynolds  and  Vonder  Haar,    1977) 


M     ms      =  AclrpclrHs  +  AcldpcldHs        cld  =  c1oud 

(2)  m,_       =  Ac]r  Nclr  ♦  Acl(j  Nc]d  c         clr  =  clear 

+  Acld(1.0-  e)Ndr 

m  -  p„i„  H.            A.,.  =  area  of  array  covered  with 

(3)  Arl,  -  ^ cui_s  cld 

cld      Hs<pcl<f   pclr> 

m  =  measured  visible  radiance 

jL  '     clr                                   p  «  .  =  bi-directional    reflectance 

«*>     Ncld=cAlrt          *Nclr  ,     „ 

cld  cloud  areas 

H.  =   solar  irradiance  reaching 


cld        clr 


s 


the  surface  (  .5  <  x  < 


(5)  Add  +  Aclp  =   1.00  .8  urn) 

in        -  N   .  "  i     3  bi-directional    reflecting 

(6)  e         =  J= -£i£  clr       of  clear  areas 


Ncld  =   sPectra1    IR  (10  <  x  < 
12  urn) 


N  «     =   spectral    IR  radiance  from 
clear  area 


m.        =  measured   IR  radiance 

e         *  IR  emissivity 

Ac^r  =  area  of  array  clear  of 
clouds 
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error  in  satellite  derived  temperature  which  is  equivalent 
to  approximately  5%  A  in  Nclr.  Ncld  will  have  a  16°C  er- 
ror if  the  Acid  has  2  %%  error  (Table  7).  This  error  can  be 
reduced  by  better  Acid  determination  using  higher  resolution 
data  and  through  use  cf  ground  truth  measurements. 

TABLE  V 

The   Errors      in   cloud   Amount    (with     A'ds   =    5%)    and   in   Cloud 
Temperature    (Nell  =    26o    K,    after   Reynolds    and   Von3er   Haar, 

1977) 


5%pclr,  10%pclr,  15%  pclr 

Uncertainty  pcld,   Ms  pcld  Ms  Pcld,   Ms 


AAcld    (Acid   =   0.26) 

0.09 

0.14 

0.20 

A Acid    (Acid  =   0.53) 

0.13 

0.21 

0.30 

AAcld    (Acid  =  0.83) 

0.17 

0.29 

0.40 

5%  ML,      '         10%  ML, 
Nclr,    z  Nclr,    £ 

24%  in  Acid        40%  in  Acid 


ANcld  16°  30° 

AZ  3   km  4.2   km 


The  bispectral  method  was  tested  with  NOAA  satellite 
coverage  in  an  area  that  had  cloud  observations  and  upper 
air  data  available.  Three  locations  were  chosen;  White  Sands 
Missile   Range    (WSMR) ,    Denver,    and  Oklahoma    City.         A    75    X    75 
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km  area  was  used  because  it  best  corresponded  to  surface  and 
upper  air  data.  Average  values  of  pcld,  pclr,  and  Hs  were 
assumed  initially  and  then  changed  as  sun  angle  and 
satellite  viewing  angles  varied  (Reynolds  and  Vonder  Haar, 
1977),  The  results  showed  a  rms  error  of  0.2  in  cloud 
amount  with  a  slight  underestimation  bias.  The  cloud  height 
error  with  cirrus  included  was  4. 6  km  with  a  -4.2  km  bias. 
Without  cirrus  cloud,  the  rms  error  improved  to  0.5  km  with 
overestimaticn   of  0.27  km. 

This  problem  with  cirrus  is  linked  to  the  use  of  0.9 
for  the  emissivity  of  all  clouds.  Because  of  the  ice  crys- 
tals in  cirrus,  its  emissivity  is  less  than  the  0.9  value. 
Two  approaches  were  considered  by  Reynolds  and  Vonder  Haar 
(1977)    to   correct  this   problem. 

One  approach  was  to  correlate  the  visible  spectral 
albedo  of  a  clcud  to  its  emissivity  through  the  blackbody 
temperature.  Fig.  20  shows  the  scatter  diagram  of  the  rela- 
tionship between  the  albedo  and  the  blackbody  temperature 
and  a  second  order  fit  which  was  used  to  determine  the  cloud 
top  temperature  for  a  specific  albedo.  Once  this  cloud  top 
temperature  is  known  the  emissivity  can  be  calculated 
using   Eq .      4,         Table   17     (derived   from   Sq.       2      of   Table    IV). 
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Fig.  21  is  the  resulting  curve  which  relates  albedo  to  eois- 
sivty.  The  curve  is  a  second  or der  equation  of  the  best  fit 
curve  measured  in  the  bispectral  method  along  with  a  slight 
adjustment  to  pclr  and  pcld.  This  cut  the  error  of  the  old 
method  in    half    (Reynolds   and   Von der  Haar,    1977) . 

The  second  procedure  for  determining  cirrus  cloud 
height  was  to  compare  channels  4  and  3  of  the  Vertical 
Temperature  Profile  Radiometer  (VTPR)  on  the  NOAA-u 
satellite.  This  led  tc  a  reduction  by  one  half  in  the  appar- 
ent error  of  the  old  method.  Table  VI  lists  the  results  of 
using  a  cirrus  correction  and  the  VTPR  approach  compared 
with   non-adjustment    for   cirrus. 

D.        PRECIPITATION 

Four  approaches  for  estimating  precipitation  from 
satellites  are  reviewed  in  this  section.  They  are  cloud 
index  (Barrett  and  Martin,  1981;  Barrett,  1981),  life  his- 
tory (Barrett  and  Martin,  1981;  Griffith  and  Hoodley,  1981), 
bi-spectral  (Liljas,  1981a,  198  1b;  Austin  and  Lovejcy, 
1981),  and  cloud  modeling  (Wylie  ,  1979;  Simpson  and  Sfiggert, 
1969) . 
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Figure  20. 


Ine  Relationship  3etween  Albedo  and  Blackbody 
naaiation.   This  was  derived  from  NCAA  satellite 
measurements  of  cirrus  clouds  ever  the  Gulf  of 
C|l^-orma  (from  Reynolds  and  Vonder  Haar, 
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Fiqure   21. 


The   Relationship   Between    Albedo   and    IB 
Emissivity   from   Empirical   Studies    (from  Reynolds 
and   Vonder  Haar,    1977). 


TABLE    VI 

9 

Comparison   of   cirrus   correction  and  VTPR    results 
Reynolds  and   Vondsr   Haar,    1977) 


(from 


CIRRUS  •CORRECTION  RESULTS 

Oate 

Uncorrected 

Si-Spectral  Cloud 

Top  Temp.    (   K) 

Cloud  Height 
Error   (Km) 

Si-Spectral   Cloud 

Top  Temp,   with 

Cirrus  Correction  ( K) 

Cloud   Heiqnt 
Error  (Km) 

Observed  Cloud  Top 
Temp.     (  K) 

S  Jan  75 

258 

♦3.2 

248 

♦1.7 

233 

8  Jan  75 

261 

♦3.3 

242 

♦  1.1 

233 

23  Oct  75 

256 

♦2.2 

229 

-  .5 

233 

9  Oct  75 

250 

RMS  Error  • 
Bias 

♦1.6 

2.6  km 
♦2.6  km 

231 

RMS  Error  -     1 
Bias            •     ♦ 

-   .4 

1     km 
48  km 

235 

VTPR  RESULTS 

Oat* 

Spot 

VTPR  Determined 
Cloud  Height  (Km) 

SR  Determined 
Cloud  Amount 

Observed  Cloud 
Heignt  (Km) 

Error  (Km) 

3  Jan  75 

33M  !<4U 

8.2 

.7 

8.7 

♦  .5 

8  Jan  75 

J3.2B   114. 3W 

8.6 

.7 

8.7 

♦  .1 

3  Jan  75 

33. 5N   115. :w 

5.7 

RMS 
Bias 

.6 

Error  •   1.75  km 
■♦1.2     km 

8.7 

♦3.0 
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1.      Cloud   Index 

Barrett  (1981)  developed  a  cloud  index  method  of 
determining  rainfall  (called  the  •Bristol  Method').  It  is  a 
manual  method  which  uses  obssrved  rain  gage  data  and  the 
cloud  index  system  (see  Section  B)  to  predict  the  amount  of 
rainfall  accumulated  over  a  spacified  time  period  (usually 
12  hours)  in  grid  squares  of  a  selected  size  (1/6°  to  1° 
grid  sguares).  Cloud  area  and  altitude  are  also  taken  into 
account  (Barrett  and  Martin,  19  81).  Fig.  22  presents  the 
flow  diagram  of  a  recent  version  of  the  Bristol  Method  with 
two    stages,      preparatory   and   operational.  Both   IS   and   VIS 

are  used  whenever  possible  to  identify  cloud  types.  Next, 
the  cloud  types  of  each  grid  square  are  assigned  an  index 
(see   Table   VII).  The    cloud    index   is      ganaraliy    "evaluated 

differently  for  ranges  of  chosen  cloud  types  and  for  differ- 
ent regions"  (Barrett  and  Martin,  1981).  Then  mean  cloud 
indices  are  calculated  along  with  mean  observed  rainfall 
amounts  for  each  cloud  field  unit.  Subsequently,  regression 
techniques  are  used  to  relate  the  observed  rainfall  to  the 
cloud  type,  and  these  help  establish  rainfall  estimates  for 
areas  of  satellite  cells  where  conventional  data  are  sparse 
(see    Fig.    23) .        Tests  show   that    the    accuracy    of  this    method 
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is    approximately  15%      for   all   cases   belonging      to  a   specific 
category    (Barrett   and   Martin,    1981). 


PflEPARATORY 
STAGES 


OPERATIONAL 
STAGES 
A 


Obtain  high  resolution 
polar-orbiting  satellite 
images  lor  lime  Itl 


Subdivide  area  into 
gacells  &  satcelle.  4 
establish  suitable 
morphoclimatic  regions 


T 


Decode  GTS 
conventional 

weather  messages 


Determine  suitable  scales  for 
cloud  type  &  cover.  &  compute 
regressions  relating  satellite 
cloud  indices  &  observed  rain 
for  each  morphoclimatic  region 


Determine  suitable  scales  to 
adjust  for  effects  of  altitude 


X 


Overlay  satellite  images  with 
grid  square  map 


Construct  simple  station  models  cepicting 
selected  weather  observations  for  all 
reporting  gacells 


Adjust  rainfall  observations  for  effects 
of  altitude 


Establish  cloud  indices  for  all  ratncloud- 
covered  grid  squares  for  time  (t| 


Calculate  mean  cloud  indices  &  observed 
rainfall  amounts  for  each  cloud  tield  unit 


Translate  satellite  cloud  indices  into  rain- 
fall estimates  via  regression  diagrams 


Compare  mapped  rainfall  field  for  time  tt) 
with  time  It  -  1 )  &  interpolate  extra 
estimates  to  allow  for  cloud  mwdvpt 


Complete  unit  time  penod  rainfall  map 


Figure  22. 


Flow   Chart   for   a   Recent   Version    of   the   Manual 
Nephanaly sis    (from  Barrstt  and    Martin,    1931). 
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TABLE   VII 

Empirically   derived    Hairfall  Probabilities   and   Intensities 
from  Satellite   Pictures    (  from    Barrett   and   Martin,    1981) 

12  3 

Assigned   probabili-  Assigned    intensities 

States  of   the  sky  ties  of   rainfall  of   rainfall    (rela- 

(nephanalysis  cloud  (relative   scale   range  tive   scale   range 

categories)  0-1.00)  0-1.00) 

Cumulonimbus  0.90  0.80 

Stratiform  0.50  0.50 

Cumuliform  0.10  0.20 

Stratocumuliform  0.10  0.01 

Cirriform  0.10  0.01 
Clear   skies 

2-      Life   History 

The  Grif fith/Woodley  (Griffith  and  Woodley,  1981; 
Griffith  et  al. ,  1978)  technique  for  estimating  connective 
precipitation  was  developed  using  a  cloud  life  history  model 
in  which  satellite  imagery  is  compared  to  radar  echo  life 
history  of  the  cloud  area  (Griffith  and  Woodley,  1981;  Bar- 
rett and  Martin,  1981).  The  basis  of  the  technique  is  the 
hypothesis  that  in  the  tropics  rain  occurs  with  cold 
(bright)  clouds  and  the  intensity  of  the  rain  is  a  function 
of  the  stage  of  development  of  cumulonimbus  clouds.  In  a 
study  of  clouds  over  Florida,  Griffith  and  Woodley  (1981) 
found  that    the      cumulonimbus   cloud  area   is      related   to    radar 
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Figure  2  3.   Regression  Diagram.  Bsed  tg  "Float"  Rainfall 

Estimates  for  Satellite  Pictures  (from  Barret- 
and  Martin,  1  981)  . 


echo  area  (see  Fig.  24).  cloud  areas  (Ac)  are  evaluated  in 
either  VIS  or  IS  satellite  imagery  and  then  compared  to  the 
radar  echo  area  (Ae)  .  Fig.  24  is  a  plor  of  the  relationship 
between  cloud  and  echo  areas  (Griffith  at  al.  ,  1978).  It  was 
empirically  derived  using  measurements  and  comparisons  of 
Florida  clouds  viewed  by  radar  and  geosynchronous 
satellites.  The  rain  cloud  threshold  (IP.)  used  in  the 
approach  was  253  K. 
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Figure   24. 


normalized  to   relative   maximum    cloud   area     (Am) 
(from   Griffith   et   al.,1978). 


This  technique  is  fully  automated  (Figure  25)  and 
has  been  tested  in  the  Florida  area  and  in  extra-tropical 
convective  areas  (Barrett  and  Martin,  1981)  .  There  are  four 
major  computer    modules  that   together    produce    a    rain   map.    The 
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programs  are  navigation,  cloud  isolation  and  tracking,  cloud 
life  history  compilation  and  volumetric  calculation  of  rain, 
and  rain  mapping  (Griffith  and  Soodley,  1981;  Barrett  and 
Martin,  1981)  . 

This  model  has  been  applied  on  several  time  scales 
and  in  different  regions  as  shown  in  Table  VIII,  The  error 
statistics  (bias  (3),  rcot  mean  square  error  (Eras} ,  and 
linear  least  squares  fit  (K,  slope,  and  intercept))  for 
these  applications  (Table  IX)  show  that  on  average  the  model 
results  agree  fairly  well  with  the  ground  truth  data 
(Griffith  and  Woodley,  1981).  The  Erms  is  small  except  for 
two  locations,  with  an  overall  average  of  0.35  millimeters 
per  hour.  The  bias  is  close  to  1.00  which  is  the  value  of 
perfect  correspondence.  In  the  linear  fits,  ail  correla- 
tions (R)  are  greater  than  or  equal  to  3.78,  where  1.0  is 
perfect  correlation. 
3.   Bispectral 

Lovejoy  and  Austin  (1979)  developed  a  technique  in 
which  visual  and  infrared  wavelength  patterns  are  used  to 
identify  clouds  with  a  high  probability  of  rain.  These 
clouds  have  large  optical  and  physical  thickness  "as 
measured   by  the   visible  albedo   in   comparison  with   their 
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Figure  25. 


The   Griff ith/Woodle?    Technique    Flow   Chart 
for   GATE     {from   Barrett   and   Sartm,    1981). 


Used 
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TABLE   VIII 

Applications  of   the   Griff it h/Woodlev   Technique    (from 
Griffith    and   Woodley,    1981) 


AREA 

Global  Atmospheric  Research  Project's 
Atlantic  Tropical  Experiment   (GATE) 

Florida  Area  Cumulus  Experiment 
(FACE) 

Equatorial  Pacific  Ocean  Climate 
Studies    (EPOCS) 

High  Plains  Experiment    (HIPLEX) 


Hurricanes 
Flash  floods 


PERIOD 
June   27   to  Sept.    20,    1974 


Selected  days  June  to  September 
1972  to  1980 

2-1/4  days  November,    1978 
30  days  November,   1979 

Selected  days.  May  to  July 
1976  to  1978 

Selected  storms  1969  to  present 

8ig  Thompson,  CO  7/31  -  8/1,    1976 
Johnstovci,   PA         July  19-20,    1977 


TABLE  IX 

The  Error  Statistics  from  the  Applications  Listed  in  Table 
VIII  (from  Griffith  and  Woodlsy,  198  1) 


LOCATION 

GATE 

HIPLEX 
(dense) 

FACE 

HIPLEX 
(sparse) 

HURRICANES 

FLASH  FLOOD 


_N  B 

53  0.84 

15  0.77 

11  0.88 

9  0.90 


1.08 
1.39 


ERM6 
0.79 
2.81 

0.42 
0.55 

0.06 
0.48 


R_  SLOPE 

0.87  0.97 

0.90  0.62 

0.95  1.18 

0.78  1.81 

1.00  0.89 

0.99  0.89 


INTERCEPT 

-1.47 

0.07 

-0.23 
-0.27 

0.18 
23.93 


PERFECT 
CORRESPONDENCE 


1.00 


0.00 


1.00        1.00 


1.00 
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height,  determined  by  the  intensity  of  the  IB  emmission" 
(Austin  and  Lovejoy,  1981).  Liljas  (1981a,  1981b)  followed 
this  with  a  simliar  type  of  bispectral  identification  of 
precipitation  area  based  on  VIS  and  IR  imagery.  This  sub- 
section will  first  review  Austin  and  Lovejoy*s  (1981) 
technigue    and   then   Liljas*^   (1981a,    1981b)     technique. 

Austin  and  Lovejoy  (1979)  first  addressed  the  ques- 
tion, "How  well  can  rain  areas  be  delineated  in  coincident 
visible  and  ir.fared  satellite  images?"  rhey  felt  that  the 
problem  was  to  find  boundaries  between  rain  and  no  rain. 
Estimation  of  rain  amount  is  assumed  to  depend  on  two  steps, 
the  delineation  of  the  rain  areas  and  estimation  of  rain 
rate  once  the  area  is  known.  The  authors  used  three  statis- 
tics to  characterize  errors  and  from  these  were  able  to 
estimate  rain  area  and  rain  rate.  The  first  step  in  rain 
estimation  was  to  map  a  frequency  matrix  for  both  rain  areas 
and  non-rain  areas  (see  Fig.  26  and  27).  The  ratio  of  rain 
to  non-rain  was  found  and  plotted  (see  Fig.  28)  along  with  a 
line  of  greater  then  50%  probability  of  rain.  From  this,  a 
rain  map  (see  Fig.  2  9)  was  produced.  The  next  step  was  to 
try   to   quantitatively      estimate    the   rain    rate.  Austin   and 

Lovejoy    (1981)       tried   several    techniques      to    do   this   but    met 
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with  little  success.  They  believe  this  was  due  to  "...the 
visible  and  IE  wavelengths  primarily  responding  to  the  rela- 
tive abundance  of  cloud  droplets  and  not  to  precipitation 
sized  particles"  (Austin  and  Lovejoy,  1981)  .  They  concluded 
that  a  good  or  even  passable  quantitative  estimate  of  rain- 
fall rate  is  not  realistic  or  practical  using  visual  and 
infrared  data.  The  authors  fait  their  study  supported  the 
hypothesis  that  visual  and  infrared  satellite  images  are 
good  for  estimating  rain  areas  but  not  for  estimating 
guantitative   rain   rates. 

Liljas'  ,(198 1a)  main  objective  was  to  apply  an 
approach  similiar  to  that  used  by  Lovejoy  and  Austin  (1979) 
for  identifying  precipitation  areas  and  intensity  using 
TISOS-N  satellite  imagery,  but  without  using  radar  to  cali- 
brate it.  He  decided  not  to  try  quantifying  intensity,  but 
to  instead  define  it  qualitative ly  by  categories  of  light  to 
heavy  rain.  This  fellows  the  work  of  ttuench  and  Keegan 
(1979)  on  the  relation  of  cloud  reflectance  and  intensity  of 
rain.  Pursuing  their  concepts,  Liljas  (1981a,  1981b)  devel- 
oped a  relative  precipitation  intensity  identification 
technique  which  relies  on  his  cloud  typing  model  (see  Sec- 
tion B,  Chapter  II)  .  It  can  be  applied  quickly  to  both 
frontal     and   convective      type      precipitation      and   gives      the 
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Figura  26. 


Frequency   Distribution  of   VIS    and   IB   Data    Pairs 
for   the   No-rain  Case     (from   Sarrett    and   Martin, 
1981). 
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Figure   27. 


Frequency   Distribution  of   18    and    VIS   Data    Pairs 
for   the   Rain    Case   at    the    Same    Location   as    Fig. 
26    (from    Earrett    and    Martin,    1981). 
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Figure  28. 


Probability  of  Rain  in  Percent  from  the  Arrays 
of  Fig.   26  and  27   with  the  50%  Boundary 
Sketched  (from  Barrett  and  Martin,  1981). 
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Figure  29. 


The    Resulting   Satellite    Pain    Sap    from    the    yse    oi 
*he   Same    Data    as    Fig.    25    and    27.       The    vertical 
lines   are  radar  rain    areas   and    horizontal    lines 
are    satellite    rain   areas    (hatched   areas    are   the 
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information  in  a  real-time  fashion-  The  first  step  is  to 
identify  the  cloud  types  (see  Fig.  9).  Next,  the  area  on  the 
graph  representing  cumulonimbus  and  niibostratus  is  sec- 
tioned (see  Fig,  10)  based  on  the  premise  that  cold,  high 
and  dense  clouds  give  the  highs st  precipitation  rate.  An 
additional  requirement  for  tha  cloud  top  temperature  to  be 
colder  than  -22°C  was  established  (Liijas,  193  1a)  and  sun- 
ported  by  investigations  of  Muanch  and  Kaagan  (1979).  Table 
X  gives  the  thresholds  determined  by  summing  the  infrared 
and  visual  digital  values  that  correspond  to  the  precipita- 
tion intensity  levels.  The  method  gave  good  information  on 
the  distribution  of  precipitation  areas  and  the  relative 
intensity  within  the  areas  for  both  frontal  and "scattered 
showers.  In  addition,  the  method  detected  "intense  precipi- 
tation in  mesoscale  cloud  clusters  between  synoptic" 
observations  (Lilias,  198  1  a) . 
U.      Cloud    Model 

Wylie  (1979)  compared  tropical  and  mid-latitude 
cloud  area  and  rain  rates  to  see  if  the  Griff ith/Wocdley 
(Griffith      et   al.,  1978)      and      other      techniques   could      be 

applied  to   non-tropical    cloud   features,    instead   of   the    trop- 
ical  cloud   systems    for  which      they   were   developed.      3ix    days 
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TABLE    X 

Threshold   Values  Describing   Precipitation    Intensity   Levels 
as  Applied   in  Fig.    10    (from   Lil jas, 198 1a) 


Chi  +  Ch4  291-310  light   rain 

311-330 
331-350 
351-370 
371-390 
>   390  very  strong  rain 


o 


f  two  different  typss  of  rain  clouds  (frontal  convective 
and  frontal  stratiform)  in  the  Montreal,  Quebec,  Canada  area 
were  chosen  for  comparison  to  the  1979  GATE  data  set  of 
Stout  et    al.     (1979)  (see    Table   711). 

The  rain  rates  of  the  two  areas  are  shown  in  Fig. 
30.  The  GATE  rain  rates  are  higher  than  "the  rain  rates  of 
Montreal,  even  at  the  colder  thresholds.  The  Montreal  rain 
rates  tend  to  vary  considerably,  especially  the  convective 
cases  (June  cases).  This  shews  that  a  scheme  based  on  the 
GATE  data  would  overestimate  rain  rate  in  the  midlatit udes, 
and  additionally,  if  one  uses  a  constant  adjustment,  there 
would   be    a   large  variance    in   the      rain      aaounts.  Table   XI 
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shows  the  comparison  for  GATE  and  Montreal  rates  of  precipi- 
tation. Wylie  (1979)  found  that  the  precipitable  water  par- 
alleled the  difference  in  the  average  rain  rate  with  some 
minor  exceptions  when  the  troposphere  was  exceptionally  sta- 
ble. He  decided  to  use  Simpson  and  Wiggert's  (1969)  one-di- 
mensional cloud  model  to  assess  the  effects  of  stability  and 
to    possibly   improve   his   results    in  the   mid-latitudes. 
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Figure  30. 


Rain   Rates  Averaged    Over   Cloud    Areas    Measured    on 
IE   Images   for    Various  Cloud   Threshold 
Temperatures     (from    Sylie,     1979). 


This  model  simulates  cloud  growth  by  using  an 
ascending  2  km  bubble.  The  model  incorporates  cloud  physics 
processes    in      the   development    of    precipitation       (Simpson   and 
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TABLE    XI 


Estimated   and  Observed   Precioit  ation   and    Stability  £diu<*t- 
ment   Factors    (after   Wylie,    1979)  * 


Prccipi- 

(Montreal/GATE) 

Predpiutkm 

Observed 

table 

Precipi-         Model 

Adjusted 

rainrate* 

water 

table          prectpt- 

satellite 

Radar 

Error 

Case 

(mm  h"') 

(mm) 

water           tat  ton 

(KTbO 

Bias 

factor 

GATE 

2-2 

54 

—                — 

.     — 

— 

— 

— 

Montreal 

29  June 

10 

38 

0-7               0-5 

24  6 

26-5 

0>» 

M 

2  June 

0-9 

22 

0-4                0-3 

8-9 

*•* 

10 

M 

26  Sept. 

0-9 

27 

0-5                0-3 

30-8 

40-3 

0-H 

1-3 

1  June 

0-4 

27 

0-5                0-3 

10-4 

7=4, 

1-4 

1-4 

20  Sept. 

0-4 

27 

0-5                00 

3-2 

f4-6 

0-2 

4-6 

16  Sept. 

0-1 

22 

0-4                00 

70 

7-2 

10 

1-0 

Average 

14-1 

17-5 

0-9 

1-7 
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Wiggert,  1969) .  Only  one  bubble  ascent  was  used  although 
normally  the  bubble  would  ascend  many  times.  As  a  result, 
the  rainfall  prediction  could  not  be  used  directly,  but  was 
nevertheless  used  as  an  adjustment  to  correct  the  Montreal 
data    for    stability. 

Table  XI  shows  the  results  of  Wylie's  (1979)  study: 
the  corrected  precipitation  of  the  satellite  closely  matches 
the  radar  measured  amount  of  rain.  Ke  noted  that  the  over- 
all  bias  was  small,  and  therefore  by  using  a  stability  cor- 
rection (calculated  from  the  one-dimensional  cloud  model), 
the  technique  to  measure  precipitation  by  satellite  in  the 
tropics   could    be   applied    to  mid- latitude   areas. 
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III.       SPADS    CLOUD    MODEL    DESCRIPTION 

A.        INTRODUCTION 

The  cloud  analysis  model  adapted  for  the  SPADS  interac- 
tive computer  is  a  combination  of  procedures  drawn  from  sev- 
eral previously  tested  techniques.  Ths  mainstay  of  the 
model  is  Liljas'  (1981a)  paral lelpipedcn  (Fig.  9  and  10) 
classification.  Added  to  this  are  a  Texture  rest  for  dis- 
criminating some  cloud  types,  a  non-linear  least  squares 
curve  fit  to  discern  cirrus  and  alto  stratus,  Harris  and  Bar- 
rett"s  (1978)  cloud  amount  calculation,  part  of  Reynolds  and 
Vonder  Haar's  (1977)  cloud  top  temperature  calculation,  and 
Liljas'  (1981a)  qualitative  precipitation  intensity  model. 
The  following  sections  describe  these  parts  of  the  SPADS 
model  and  explain  why  they  were  selected.  The  last  section 
of  the  chapter  summarizes  with  the  flow  chart  of  The  model's 
computer   program. 

3.       CLOUD    TYPING 

The  approach  selected  for  cloud  Typing  was  developed  and 
tested  by  Liljas  (1981a).  It  is  a  mult ispectral  method 
that   uses      three  channels   to     type   clouds      and    differentiate 
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between   land     and   water-  Ths    SPADS      cloud   model      uses   two 

channels,  thereby  becoming  bispectral  (the  near-IR  channel 
was  not  used  due  to  non-availability  in  GOES  VISSR  data). 
Fig-  10  depicts  the  bispectral  model  iicorporat ed  in  the 
SPADS  Cloud  Model.  Both  the  visual  and  infrared  thresholds 
shown  in  Fig-  10  were  converted  from  riros-N  Advanced  Very 
High  Resolution  Radiometer  (AVHRR)  to  GOES  VISSR  data 
counts.  Liljas*  visual  data  brightness  counts  ranged  from  0 
to  255  while  the  SPADS  counts  range  from  D  to  63,  A  four  to 
one  linear  relationship  exists  in  the  VIS  channel  and  there- 
fore a  direct  conversion  of  thresholds  was  performed.  The 
IR  thresholds  were  converted  to  temperatures  (K)  using  an 
August  1980  TIROS-N  AVHRR  tei perature  calibration  table. 
This  gave  a  rain  cloud  threshold  of  about  251  K  which  is  in 
agreement  with  the  G rif f ith/Woo dley  (Griffith  and  Woodiey, 
1981)  and  the  Liljas  (1981a)  rain  cloud  thresholds.  These 
temperature  thresholds  were  converted  to  GOES  IR  brightness 
counts.  The  transformed  thresholds  used  in  the  SPADS  cloud 
model   are    depicted    in   Fig.    31. 
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One  of  the  reasons  for  selecting  the  Liljas  model  is 
that  it  was  tested  or  an  interactive  computer  siiailiar  to 
SPADS.  In   addition,        the   test     area     selected   by      Liljas 

(1981a)  included  land  and  ocean  areas  along  with  many  types 
of  clouds  associated  with  mesoscale  and  synoptic  scale  fea- 
tures. A  complete  description  with  color  pictures  is  in 
Liljas' (  1981a)  technical  report.  In  addition,  a  combined 
approach  using  texture  and  bispectral  counts  was  suggested 
by  Liljas  (1981a),  Harris  and  Barrett  (1978)  and  Reynolds 
and  Vender  Haar  (1977).  The  U.S.  Air  Force1 s  3DNEPH  program 
used  both  these  approaches,  but  documentation  of  thresholds 
and  tests  was  not  available  to  us.  Thus,  the  3DNEPH  program 
could   only    be   used   as    a    qualitative   guide. 

The  GOES  IR  resolution  is  approximately  4  n  mi.  at  the 
subsatellite  point,  but  for  the  visual  channel  a  higher  res- 
olution of  0.5  n  mi.  is  available.  This  VIS  resolution 
affords  an  opportunity  to  get  a  better  estimate  of  cloud 
cover.  All  the  studies  reviewed  highly  recommended  the  use 
of  an  array  cf  bright ress  counts  to  give  an  improved  repre- 
sentation of  the  cloud  picture.  The  0.5  n  mi.  resolution 
allowed  an  8X8  array  of  visual  pixels  per  IR  pixel.  For 
each  grid  space  an  average  visual  brightness  count  is 
calculated    (Eg.    1    ,    Table    II). 
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Figure   3  1. 


Tvo    Dimensional  Clou  d .  Tyomg   iraoh    Using   GOES    IB 
and    VIS    Satellite    Digital  Data    (after    Llljas, 
1981a)  . 


At  the  same  time  as  the  visual  average  brightness  is 
baing  calculated,  the  standard  deviation  of  the  cloudy  vis- 
ual counts  per  grid  space  is  datermined.  The  standard  devi- 
ation (Eg.  4,  Table  I)  is  the  statistic  used  to  represent 
texture.  As  suggested  and  tested  by  Harris  and  Barrett 
(1978),  the  standard  deviation  is  used  to  discriminate 
between  stratiform  and  cumulus  humulus,  small  and  large 
cumulus  congestus,  and  altostratus  and  cirrus.  The  initial 
standard  deviation  threshold  values  for  this  discriminant 
analysis  were  approximated  by  using  information  in  Fig.  6 
and  8.  An  additional  test  was  used  first  for  the  cirrus/al- 
tostratus  decision.  A  second  order  polynomial  least  squares 
fit  was  applied  to  the  cirrus  (c)  points  and  also  to  the 
altostratus  (d)  points  in  Fig.  10.  Tha  calculated  curves 
are  depicted  in  Fig.  31,  with  the  associated  equations  in 
Table  XII.  Entering  the  equations  with  the  average  visual 
brightness,  the  resulting  IB  value  is  compared  to  the  actual 
IR  to  determine  cloud  type.  If  the  calculated  value  falls 
between  the  two  curves  then  an  additional  test,  using  the 
standard  deviation  threshold,  is  used  to  discriminate  cloud 
type. 


82 


C.       CLOUD    TOP    TEMPERATURES 

The  Reynolds  and  Vonder  Haar  (1977)  •  study  is  used  as  a 
basis     for        estimating      cloud      top        temperatures.  Their 

calculation  is  based  en  two  equations  (23.  1  and  2,  Table 
IV) ,  one  for  cloud  amount  and  the  other  for  the  cloud  top 
temperature.  The  resolution  of  the  data  used  was  4  n  mi. 
but  the  authors'  analysis  suggested  an  improved  resolution 
would  give  better  results.  As  mentioned  above,  this  led  to 
the    choice   of   0.5   n   mi.    visual  data   in   this    study. 

The  first  step  in  determining  the  cloud  top  temperature 
is  to  calculate  the  a  noun t  of  cloud  per  IB  pixel.  The  cloud 
decision  was  based  on  comparing  the  visual  digital  count  of 
each  pixel  in  the  8X8  grid  to  a  no  cloud  threshold.  This 
threshold  value  (initially  5  VI S  counts)  was  derived  from 
Fig.  31  (after  Liljas,  1981a)  where  it  is  the  lowest  VIS 
count  to  have  clouds.  A  cloud  threshold  test  was  also  used 
by    Harris      and   Barrett    (1978).  Their   threshold      value   was 

derived  by  comparing  the  visual  image  with  a  computer  pro- 
duced density  slice,  but  neither  the  value  nor  the  formula 
to  obtain  the  value  were  included  in  their  report.  The  U.S. 
Air  Force's  3DNEPH  program  based  the  cloud  decision  on  the 
amount      of    visual      grayshade    difference      between   a      3    n      mi. 
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TABLE   XII 

SPADS   Cloud  Model's   Cloud   Types,    Least   Squares. Fit 
Equations,    and    Standard   Deviation   Values    (SIGi) 


Cirrus    (1)*  -  Cs 

a.  y  >   158.094  -  6.6076  x  +  4.159  x2  ** 

b.  a  <_  SIG3 
Altostratus  (2)  -  As 

a.  y  £  115.24  +  .5246  x  +  .0395  x2  ** 

b.  a  >  SIG3 

*  The  number  in  the  parens  is  the  computer  cloud  type  identifier, 

^r  y  =  IR  digital  value 

x  =  average  Visual  digital  value 

Fog/stratus  (3)  -  ST 

Thick  Fog/stratocumulus  (4)  -  Sc 

a  <_  SIG  1 
Cumulus  humilis  (5)  -  CuHu 

a  >  SIG  1 
Small  cumulus  congestus  (6)  -  CuCong 

a  £  SIG  2 
Large  cumulus  congestus  (7)  -  CuCong 

a  >  SIG  2 
Nimbostratus/multilayered  (8)  -  Ns 
Cumulonimbus  (9)  -  Cb 


an 


pixel  and  a  background  brightness  (Fig.  15,  refer  to  Chap- 
ter II  for  more  complete  details).  This  program  uses  an 
extensive  global  background  brightness  file  that  is  updated 
regularly  (Fye,  1978).  The  development  and  maintenance  of 
such   a  file   for    SPADS   is    not   feasible   at   this   time. 

The  average  cloud  amount  is  determined  by  summing  the 
number  of  cloud  decisions  in  the  grid  space  and  dividing  by 
the  total  number  of  pixels  per  grid  space.  This  calculation 
followed  the  Harris  and  Barrett  (1  978)  Eg.  3  (Table  I)  . 
This  cloud  amount  fraction  is  then  used  in  the  Reynolds  and 
Vender  Haar    (1977)       eguation   for    cloud   top    IR,  radiance    (   Eq. 

9 

2,  Table  IV).  On?  of  the  reasons  this  cloud  amount 
calculation  is  used  instead  of  Reynolds  and  Vonder  Haar's 
eguation  was  that,  in  the  latter,  three  additional  constants 
would  have  to  be  calculated  and  verified,  which  is  beyond 
the  scope  of  this  thesis.  The  authors  felt  that  it  was 
essential  in  using  their  cloud  amount  equation  that  these 
constants  be  well  determined  (Reynolds  and  Vonder  Haar, 
1977) . 

The  next  parameter  of  the  cloud  top  radiance  equation  to 
be  determined  is  the  clear  area  IB  spectral  radiance.  Fleet 
Numerical    Oceanography      Center    ( FNOC)         analyses   of      surface 
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temperatures      (K)  are      obtained        and        a      representative 

temperature  profile  matching  the  location  of  the  IR  image 
center  point  is  selected.  These  temperatures  are  converted 
to  IR  radiances  via  Planck's  function  and  used  in  the  cloud 
top    temperature   equation, 

Emissivity  of  the  clouds  also  must  be  determined. 
Reynolds  and  Vonder  Haar  (1977),  among  others,  found  that 
the  emissivity  is  nearly  constant  (approximately  equal  to 
0.9)  for  all  clouds  except  cirrus,  for  which  it  is  approxi- 
mately 0.55.  In  the  program,  cirrus  clouds  are  identified  by 
the  cloud  typing  and  standard  deviation  nodules  so  that  the 
lower  emissivity  can  be  used  in  the  cloud  top  temperature 
equation.  In  addition,  an  emissivity  of  1.0  was  used  for 
cumulonimbus  and  nimbostratus  based  on  the  assumption  that 
these   clouds    do   not   allow    any   surface   radiance    through. 

The  cloud  top  radiance  (Ncld)  is  calculated  using 
Reynolds  and  Vonder  Haar's  (1977)  equation  (Zq.  2,  Table 
III)  which  consists  of  the  surface  temperature  field 
obtained  from  ?NOC  (converted  to  an  IR  radiance  via  Planck's 
law),  the  emissivity  (cirrus  or  non-cirrus),  the  satellite 
IR  radiance  (Ml)  value,  and  the  fractional  amount  of  cloud 
(Acid)  .    Once   Ncld   is    calculated,    a   conversion    to    a   cloud   top 
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temperature  is  made  via  the  inverse  of  Planck's  function. 
This   temperature     can   then     be  compared     to   the      latest   FNOC 

« 

upper  air  data  analyses  to  find  the  corresponding  height  in 
millibars  (rab)  of  the  cloud  top.  The  upper  air  data  are 
available  at  ten  standard  levels  consisting  of  1000,  900, 
850,  700,  500,  400,  300,  250,  150,  and  100  rab. 

D.   PRECIPITATION  INTENSITY 

The  precipitation  intensity  identification  technique 
selected  is  the  Liljas'  (1981a)  method.  It  is  an  extension 
of  his  cloud  typing  i^thod  (sea  Chapter  II)  and  therefore 
easily  fits  into  the  overall  SPADS  cloud  model  without,  add- 
ing much  time  to  the  process.  Sylie's  (1979)  cloud  model,  as 
discussed  in  Chapter  II,  used  the  Griff ith/Wocdley  (Griffith 
§1  §.!•*  1973)  technique  with  a  stability  correction  which 
gave  good  indication  cf  rain  rates.  One  problem  with  it  is 
that  the  Griff  it h/ Woo dley  technique  is  based  on  the  develop- 
ment and  identification  cf  convective  cloud  systems*  life 
history  through  use  cf  radar  echoes  not  available  to  this 
study.  Another  problem  with  using  this  scheme  is  the 
computer  run  time  it  would  add  to  the  SPADS  cloud  model. 

Although  time  is  an  important  factor,  it  is  not  the  only 
reason  for  selecting  the  Liljas  method.  Liljas  (1981a) 
reported  good   success  with  this  technique,    obtaining  good 
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indications  of  relative  precipitation  over  land  and  sea 
areas  as  verified  by  synoptic  reports.  He  felt  that  this 
technique  would  be  useful  for  sea  areas  where  no  synoptic 
observations  were  available.  In  addition,  relative  precipi- 
tation from  mesoscale  cloud  clusters  was  detected  between 
synoptic  observations.  This  identification  scheme  was  modi- 
fied for  use  in  the  SPADS  cloud  model  by  red-icing  the  number 
of  precipitation  intensity  categories  froa  six  (see  Fig.  9) 
to  three.  Fig.  32  depicts  the  new  categories  of  precipita- 
tion types  representir.g  light,  moderate,  and  heavy  rainfall 
and  Table  XIII  includes  the  equation  for  the  lines  separat- 
ing the  categories.  The  light  rainfall  category  includes 
the  requirement  that  cloud  top  temperatures  be  less  than 
251     K. 

E.        FLOW    CHART 

The  SPADS  Cloud  Model  (Fig.  33)  begins  with  obtaining 
the  initial  infrared  and  visual  data  fields  from  the  GOES 
satellite  and  conventional  surface  and  upper  air  data  from 
FNCC.  The  IB  and  VIS  satellite  data  are  collocated  both  in 
space  and  time.  The  upper  air  and  surfaoe  temperatures  are 
selected  to  correspond  to  the  I, J  position  closest  to  the 
center  point  of   the  image  from  the  latest   FNOC  analysis  on 
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the  same  date  as  the  satellite  pictures.  Next,  the  SPADS 
cloud   model   computer    program    (see   Appendix    A)    is   run. 

As  shown  in  Fig.  33,  the  first  step  of  the  program  is  to 
calculate  three  values  based  on  the  8X8  grid  of  VIS  pixels; 
the  average  VIS  brightness,  standard  deviation,  and  cloud 
amount  as  discussed  in  Section  B  of  this  chapter.  Once  the 
values  are  obtained,  the  cloud  typing  module  {based  on 
Liljas,  1981a)  is  called.  Two  tests  are  performed  within 
this  module;  on*  is  based  on  comparison  of  the  IR  and  aver- 
age VIS  counts  and  the  other  supplements  the  first  by  using 
texture  (standard  deviation)  to  complete  the  cloud  typing 
task.  If  nimbostratus  or  cumulonimbus  clouds  are  identi- 
fied, the  precipitation  module  is  called  and  it  determines 
the    intensity    (gualitatively)    of    the    precipitation. 

The  cloud  top  temperature  module,  in  which  the 
temperature  and  the  height  of  the  clouds  are  calculated,  is 
called  for  all  cloud  cases.  The  FNOC  surface  temperature 
value  is  used  in  the  cloud  top  temperature  calculation, 
while  the  FNOC  upper  air  data  are  used  in  the  height 
calculation. 
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Figure  32. 


Graph  of.SPADS   cloud    Model   Precipitation 

Intensities,  fc-wi^iauw 
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TABL2  XIII 
SPADS  Cloud  Model  Precipitation  Intensity  Categories 


Equation 


suml  =   IR  radiance  +  average  VIS   brightness 


No   rain   (<£) 


suml  <    184 


Light   rain    (1) 


184   <    suml    <    195 


Moderate   rain   (2)  195   <   suml  <   224 


Heavy  rain    (3) 


224   <   suml 


The  final  step  is  processing  the  results.  Cloud  types 
(Table  XII),  precipitation  intensities  (Table  XIII),  cloud 
top  temperatures  (K)  ,  and  cloud  top  heights  (mb)  are  stored 
in  a  permanent  file  f cr  future  display.  The  average  visual 
brightness,  standard  deviation,  and  the  amount  of  cloud  that 
correspond  to  each  IR  pixel  are  printed  for  reference  and 
validation    purposes. 
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INPUT 
GOES  IB  Satellite  Image 
GOES  VIS  Satellite  Image 
FNOC  Temp  (SFC) 
FNOC  Temp  (UA  mandatory) 

•" 

CALCULATE  (8x8  sq.  VIS) 

1.  Average  VIS  Brightness 

2.  Standard  Deviation 

3.  Cloud  Amount 

i 

* 

1 
p. 

DETERMINE  CLOUD  TYPE 
IE  &  AVG  VIS  -Thresholds 
TEXTUPE  (standard  deviation) 

> 

r 

> 

t 

Bispectral  Calculation 

of  Cloud  Top,  Temp 

and  Height 

Determine  Precip 
Intensity 
(Ns,  Cb) 

i 

* 

Results 

Figure  33.   The  SPADS  Cloud  Modal  Generalized  Flow  Chart. 
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IV.       THE   MODEL    TEST    AN D    RESULTS 

A.  INTRODUCTION 

A  test  was  conducted  to  determine  if  the  SPADS  Cloud 
Model  concepts  were  valid.  The  test  was  designed:  (1)  to 
capture  a  small  area  cf  coincidant  IR  and  Vis  satellite  data 
in  the  mid-latitudes  and  to  obtain  FNOC  fields  for  the  same 
date-time  group  and  location  as  the  satellite  data,  (2)  to 
run  the  SPADS  Cloud  Model  using  these  data,  and  (3)  to  ana- 
lyze the  results  of  the  model  to  determine  the  validity  of 
the  model's  analyses.  The  following  sections  describe  the 
input  to  this  test,  test  procedures,  results,  and  subsequent 
iterations  and  adjustments  made  to  improve  the  model's 
performance. 

B.  DATA 

GOES  IR  and  VIS  data  were  captured  and  archived  on 
November  9,  1982  at  1315  local  time,  far  a  site  centered  at 
location  latitude  44° N,  longitude  14  1°W.  The  corresponding 
FNOC  data  fields  were  also  saved  for  the  same  date  from  1200 
GMT  data.  Ficj.  34  is  the  GOES  IR  image  from  which  an  area  of 
approximately  460  X  460  n  mi.  was  selected  for  analysis. 
This   afforded    4096   IR   pixels    for    evaluation.    Fig.      35    is   the 
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IB  satellite  picture  of  the  selected  area  with  the  center 
point  matching  that  in  Fig-  34.  The  IR  data  resolution  was 
approximately  7  n  mi.  at  this  latitude  and  4  n  mi.  at  the 
subsateilite  point.  The  visual  data  resolution  used  was  0.9 
n  mi.  (0.5  n  mi.  at  subsateilite  point)  .  The  visual  area 
(see  Fig.  36)  is  the  same  area  as  the  infrared  giving  an  8X3 
visual  grid  of    pixels    (64    points)    par   IS   pixel. 

The  FNOC  vertical  temperature  profile  was  obtained  to 
match  the  center  point  of  the  satellite  pictures  and  the 
corresponding  heights  (mb)  are  listed  in  Table  XIV.  Only  one 
set  of  temperatures  was  used  because  the  resolution  of  the 
FNOC  fields  (63  X  6  3  field  with  38  1  ki  mesh  length)  is 
approximately  the  same  as  the  area  of  the  satellite 
pictures    used    in   this   test. 

C.        TEST    PROCEDURES 

The  following  list  is  the  set  of  procedures  for  this 
test   of   the   SPADS  Cloud    Model.    The  steps   taken    were   to: 

1.  Capture   coincident    VIS   and    IR  satellite    images. 

2.  Receive    surface      and    upper      air    temperatures      from 
FNOC. 

3.  Manually      analyze   images      for   cloud      boundaries    and 
types. 

4.  Select   the   temperature   profile   closest    to    satellite 
image. 
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5.  Run    SPADS   Cloud    Model   computer   program. 

6.  Output   results   in  a    format    for   analysis. 

7.  Compare  manual  tc  computer    output. 

8.  Adjust  computer    model  as   appropriate. 

9.  Repeat    steps   5,     6,    and  7. 

The  manual  nephanalysis  was  completed  by  Dr.  C.H.  Washr 
Assistant  Professor  in  Meteorology  at  the  Naval  Postgraduate 
School.  He  used  both  the  IR  and  Visual  images  (Fig.  34,  35, 
and  36)  to  discern  the  cloud  boundaries  and  types  seen  in 
Fig.    37    (see  Table  XV   for    symbol    definition) 

During  the  first  computer  run,  a  cloud  threshold  of  5 
visual  digital  counts  was  used.  This  value  corresponded  to 
the  lowest  visual  threshold  that  could  have  clouds  (see  Fig. 
33)  .  After  analysis  of  the  first  run,  however,  a  threshold 
of  20  counts  was  selected.  This  adjustment  accounts  for  our 
uncertainty  in  the  GCSS  VTSSR  visible  channel  calibration, 
as  compared  to  that  for  the  TIROS  AV HRR  visible  channel.  In 
normal  use,  VISSR  is  a  relatively  gualitative  radiometer, 
and  the  engineering  data  needed  to  convert  values  to  high- 
guality  calibrated  radiances  were  not  readily  available. 
Two  sets  of  standard  deviations  were  used  and  are  presented 
in    Table    XVI. 
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Figure  3U. 


Infrared    Satellite   Image,    Center    Point   - s   at 
Latitude    44°N,    Longitude    141°w.  ~"      ~ 
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Figure  35. 


Infrared    Satellite   Image    Selected    from   F 
with  an   Area    of  460    X  460    n    mi.     (center 
at   latitude    44°N,    longitude    141°W). 
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TABL2   XIV 

FNOC   TemDerature   Profile    for    12  00  SMT,    November   9,    1982   at 
latitude   445Nr    longitude    1U1°W 


Standard  level   (mb)         Temperature  (°K) 


100 

212. 

150 

214. 

200 

216.8 

250 

221.8 

300 

229.7 

400 

244.6 

500   • 

254.3 

700 

267.7 

850 

275.5 

1000 

283.4 

1010  (SFC) 

289.7 

99 


Mid-Cloud 

AS 
Possible  Ac 


Figure  37. 


Manual  Clcud    Analysis  of    Fig.     35    and    36    (see 

Table   XV    for    the   symbol   definitions). 
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TABLE    XV 
Definitions  of   Symbols   Used   in   Fig.    37 


Nimbostratus 

-  Multi 

pie  1 

ayers 

NS-ML 

Al tcstratus 

AS 

Altocumulus 

AC 

Cirrus 

Ci 

Cumul  us 

Cu 

Stratus 

St 

Stratocumulus 

SC 

Broken 

TABLE 

XVI 

0 

Standard   Deviation    (SIGiJ    Values   Used   for    Run    1    and   2    (see 

Table   XII    tor   how    they    are   used) 


SIG    (i) Test    1 Test    2 

SIG   1  4  2.5 

SIG  2  20  20 

SIG   3  3  1.5 


D.        RESULTS 

The  following  subsections  discuss  the  results  of  the  two 
computer  runs  as  compared  to  the  manual  analysis.  The 
results    are   discussed   by   the   type   of    product    and   run    number; 
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where  the      products  are      cloud   type,         precipitation,      cloud 
amount,        cloud   top      temperature      and      height,      and     texture 
(standard   deviation)  . 
1#      Run  J. 

a-      Cloud  Typing 

The  overall  general  patterns  of  the  cloud  types 
(see  Fig.  38)  compared  well  to  the  manual  analysis  (see  Fig. 
37).  The  observed  cloud  structure  is  most  easily  discussed 
by  considering  the  images  as  four  quadrants;  northeast  (NE) , 
northwest  (NW)  ,  southwest  (SW)  ,  and  southeast  (SE)  with 
north   being   the   Top   of  the   image. 

The  NE  quadrant  has  mid-level  clouds  such  as 
aitostratus  (see  Fig.  37).  The  SPADS  cloud  model  did  well 
here  indicating  predo mina tly  mid-level  clouds  (see  Fig.  38). 
The  model  also  matched  the  manual  analysis  in  the  NW  quad- 
rant. This  area  has  high  cloud  tops  of  nimbostratus  and 
embedded   cumulonimbus. 

The  SW  quadrant  has  a  mixture  of  mid- level  and 
multi- layered  thick  clouds  and  again  the  cloud  model  is  in 
agreement  with  the  manual  analysis  (see  Fig.  37  and  38).  The 
objective  analysis  shows  more  detail  than  the  manual 
analysis    because      the  computer  program      is    working      pixel    by 


102 


pixel   while  the   manual     analysis      cm   only    resolve   groups   of 
pixel   defining   cloud    areas. 

In  the  SE  quadrant,  the  SPADS  cloud  model  could 
not  resolve  the  difference  between  low-lsvel  (stratus)  and 
mid- level  clouds.  The  boundary  between  these  is  located  on 
the  upper  border  of  the  SE  quadrant.  Another  discrimination 
problem  is  batween  stratus  and  cumulus  in  the  lower  right 
quadrant.  The  amount  of  cumulus  humulus  that  seems  apparent 
in  the  visual  satellite  picture  (see  Fig.  36)  and  the  manual 
analysis  (see  Fig.  37)  did  not  appear  in  the  computer 
analysis  (see  Fig.  38).  This  problem  is  related  to  the  tex- 
ture threshold  tests,  as  both  stratus  and  cumulus  humulus 
have  the  same  cloud  type  thresholds. 
b.      Precipitation 

The  patterns  of  rain  intensity  (see  Fig.  39)  in 
the  NW  and  SW  quadrants  are  based  on  the  premise  that  the 
thicker  and  colder,  therefore  brighter,  clouds  (compare  Fig. 
39  to  Fig.  35  and  36)  have  a  greater  probability  of  rain  and 
a  higher  intensity  of  rain.  Along  with  this,  is  the  knowl- 
edge that  the  approximate  threshold  for  rain  clouds  is 
colder  than  -22  °C  (Liljas,  T981a;  3arrett  and  Martin, 
1981).      3oth   computer   runs   gave      the    same    agreement    with   the 
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Figure  33.   Cloud  Types  as  Analyzed  by  tha  SPADS  Cloud  Model 
in  Run  1  (see  Tabic  XII  for  numeric 
definitions)  . 
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manual  analysis  because  the  areas  of  rain  clouds  did  not 
change.  At  this  time,  no  surface  data  are  available  for 
ground-truth   verification  of  the   rain    patterns. 

c.  Cloud   Amount 

The  cloud  amount  for  the  first  run  was  100%  for 
the  entire  analysis  area.  The  visual,  18,  and  manual  pic- 
tures all  shewed  clear  and  scattered  cloud  areas  in  the  SS 
guadrant.  In  the  manual  analysis  (see  Fig.  37),  this  quad- 
rant has  these  sections  marked  Cu,  St/So,  .  and  (£7)  As  (see 
Table  XV  for  symbol  definitions!  .  Further  research  into  the 
cloud  threshold  for  GOES  satellite  imagery,  combined  with 
analysis  of  the  high  resolution  visual  data  used  in  this 
test,  suggested  a  better  cloud  threshold  value  would  be  20 
visual  counxs  (Muench  and  Keegan,  1979,  and  see  discussion 
in    Section   C) .  The  results   from   using      this   threshold   are 

discussed   in  the   Run    2  section. 

d.  Cloud   Top   Temperatures  and   Heights 

The  cloud  top  temperatures  (see  Fig.  40)  and 
heights  (see  Fig.  41)  look  reasonable  for  the  clcud  -ypes 
identified.  The  NE  quadrant  had  temperatures  ranging  from 
260  to  250  X  and  a  height  of  a pproixmately  500  mb,  these 
matching      the    type      of  mid-level      cloud.  The    NW      quadrant 
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Figure   39. 


Qualitative   Precipitation   Intensity    from   Run    1 
of   the    SPADS    Cloud   :iodel     (see    Table    XIII    for 
numeric   definitions)  . 
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clouds  are  mostly  nim tostratus  with  high-level,  cold  cloud 
tops.  The      height        is      approximately      3C0      mb        and      the 

temperature  is  230  K.  Again,  agreement  is  good.  The  mixed 
clouds  in  the  SW  guadrant  also  have  mixed  height  (mid-level 
and  high-level)  and  temperatures  which  compare  well  (see 
Fig.    40    and   41) 

The  SE  guadrant' s  heights,  ranging  from  850  to 
700  mb,  define  the  area  reasonably  well  (see  Fig.  41).  The 
division  for  the  low  and  mid-level  clouds  (see  Fig.  37) 
appears  as  the  700  mb  contour  line  (see  Fig.  40).  The  cloud 
temperature  field  (see  Fig.  40)  defines  the  division  of  low 
and  mid-level  clouds  less  clearly,  but  partly  cloudy  areas 
(275  K  contour) are  indicated.  The  cirrus  and  aitostratus 
clouds  are  too  cold  and  high,  as  compared  to  the  FNOC 
temperature  profile.  This  error  appears  to  be  in  the  appli- 
cation of  the  bispectral  radiance  eguation  (Eg.  4,  Table  IV) 
from  Reynolds  and  Vender  Haar  (1979)  in  partly  cloudy 
regions.  The  coefficients  in  that  eguation  depend  stror.gly 
on  the  radiance  calibration  of  the  radiometer  used.  The 
assumed  calibration  ccefficeints  used  in  this  study  are 
clearly  inadeguate  and  must  be  improved  through  empirical 
validation  of  GOES  VISSH  calibration  using  ground  truth 
data. 
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e.      Clcud   Texture    (Standard  Deviation) 

The        cloud  texture         (standard  deviation) 

calculation  results  are  presented  in  Fig.  42.  Overall,  they 
show  an  expected  pattern;  mixed  and  broken  cloud  areas  (in 
the  Sfl  and  SE  quadrants)  have  the  highest  variation  and  uni- 
form stratus  type  cloud  cover  areas  (in  the  NW  and  HE  quad- 
rants) have  the  smallest  variation  (compare  Fig.  36  and  42)  . 
The  highest  deviations  (approximately  '4.0)  were  located  in 
the  SE  guadrant  in  the  areas  discussed  previously  under 
cloud  amount.  After  contouring,  there  appeared  to  be  a  defi- 
nite boundary  between  the  highly  variable  region  (SE  quad- 
rant) and  the  nearly  uniform  areas  (NE  quadrant,  see  Fig. 
U2) .  The  SE  quadrant  had  been  manually  typed  as  cumulus, 
with  broken  areas,  and  stratus/ stratocumulus  clouds  (see 
Fig.  37)  .  After  evaluation,  the  standard  deviation  threshold 
value  for  discerning  cumulus  humulus  was  adjusted  along  with 
the  value  for  altostratus.  There  did  not  appear  to  be  obvi- 
ous patterns  of  different  size  cumulus  congestus  and  there- 
fore -his  standard  deviation  threshold  was  net.  adjusted. 
The  standard  deviation  definitely  shows  promise  in 
discriminating    between  cumuloform   and   stratus    cloud   types. 
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2.   Run  2 

The  results  from  the  second  run  are  presented  below 
with  the  revised  cloud  and  standard  deviation  thresholds. 
Only  regions  where  changes  from  Run  1  occurred  are  discussed 
in  detail. 

a.  Cloud  Typing 

The  cloud  typing  re suits  wers  changed  somewhat 
because  of  the  new  cloud  threshold  and  standard  deviation 
thresholds  (see  Fig.  43)  especially  in  the  lower  right  quad- 
rant. Here  cumulus  humulus  and  altostatus  types  were  deli- 
neated and  the  patterns  appear  closer  to  those  in  the  manual 
analysis. 

b.  Cloud  Amount 

The  cloud  amount  patterns  (see  Fig.  44)  are  more 
in  line  with  the  manual  analysis  and  what  can  clearly  be 
seen  on  the  satellite  images.  It  appears  thatthe  cloud 
threshold  of  20  is  better,  but  this  may  require  additional 
study. 

c.  Cloud  Top  Temperature  and  Height 

The  problem  in  obtaining  The  correct  cloud  top 
temperature  and  height  in  a  partly  cloudy  situation  using  a 
bispectrai  technique  still  remained.  This  problem  stems 
from  uncertainty  in  calibration  of  the  GOES  VISSR.  The  VISSR 
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Figure  u  3. 


Cloud   Types    as   Analyzed   by   tha    SPADS    Cloud 
Model,    Run  2     {see   Table   XII    for    numeric 


definitions)  . 
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infrared  channel  is  not  actively  calibrated  and  extensive 
ground  truth  is  required  to  accurately  interpret  its  data  in 
terms  of  temperatures  (Maul,  198  1).  In  the  "present  program, 
temperature  calibration  would  enter  into  the  coefficients  of 
the  cloud  height  3q.  4  (Table  17)  and  directly  into  the 
cloud  top  temperatures.  Time  was  not  available  in  this  study 
to  make  this  adjustment. 
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V.       SUMMARY 

A.        DISCUSSICN 

This    chapter   summarizes  the    thesis   accomplishments,      the 
problems    encountered,    and   suggestions   for   future   research. 

1 .      Accomplishments 

The  first  accomplishment  was  to  complete  an  overview 
of  cloud  and  precipitation  studies  with  an  emphasis  given  to 
those  using  an  interactive  computer  system.  Next,  useful 
techniques  and  ideas  were  selected  from  these  studies  and 
used  to  design  the  SPADS  cloud  and  precipitation  model.  This 
SPADS  Cloud  Model  was  then  coded  and  established  on  the 
SPADS  system.  Finally,  two  test  runs  on  a  test  case  were 
completed    with   good   success. 

Within  an  hour,  detailed  information  was  produced  by 
the  SPADS  Cloud  Model  through  concurrent  analysis  of  high 
resolution  visual  and  infrared  iigital  satellite  data.  Over 
a  maritime  region,  the  SPADS  Cloud  Model  estimated  various 
cloud  types,  precipitation  intensity,  cloud  amounts,  and 
cloud   top   temperatures  and   heights. 
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2 .      Problems 

Several  problems  were  encountered.  The  cloud  amount 
threshold  needed  adjustment  and  two  tests'  of  different 
thresholds      were      accomplished.  The      threshold   of    20      was 

satisfactory   but   further    verification   is   suggested. 

Another  area  cf  concern  is  the  threshold  values  for 
the  texture  zests.  Two  thresholds  were  adjusted  in  the  tests 
bur  all  three  need  further  refinement.  Ihe  separation  test 
(the  least  square  fit  equations)  for  cirrus  and  altostratus 
could  not  be  substantiated  in  this  data  set  and  will  also 
need   verification   in    future  runs. 

An  additional  problem  is  in  the  cloud  top 
temperature   and   height  calculations.  As    discussed    before, 

the  temperatures  are  roc  cold  and  the  heights  too  high  in 
certain    partly      cloudy  areas.  The    bispectral      equation   to 

obtain  these    values   needs   further  study. 

All  three  of  these  problems  are  related,  at  least  in 
part,  to  the  problem  of  inter-calibration  of  AVHRR  and  GOES 
VISSE   visible   and   infrared  data. 
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3 .      Recomendatio n  s 

The  following  are  reco mmendations  and  suggestions 
for    additional  research.    The   SPADS   Cloud   Model   should    be: 

a.  converted  to  use  temperatures  for  infrared  thresholds 
and  albedoes  for  visual  thresholds  to  allow  easier  con- 
version between  different  satellites  (polar  and  geosta- 
tionary) and  seasons  (this  will  require  development  and 
addition  of  sensor  calibration  tables  and  modules  to 
the    SPADS   Cloud    Model). 

b.  converted  tc  a  variable  resolution  so  that  larger  areas 
could  be  analyzed  in  the  same  time  frame  but  also  keep 
the    high   resolution   capability. 

c.  tested  and  verified  in  various  regions  by  using  surface 
observations,  rawinsondes,  aircraft  platform 
observations,  ship  observations,  radar,  and, if  availa- 
ble,   the   U.S.    Air  Force   3DNEPH  cloud   analysis. 

Digital  satellite  data  are  currently  underutilized, 
except  perhaps  for  sea  surface  temperature  and  temperature 
sounding  applications.  Most  visual  and  infrared  satellite 
data  are  treated  as  images  and  are  still  analyzed  manually. 
Now,  with  the  development  of  SPADS  and  similiar  systems,  new 
satellite  products  can  be  produced  in  regional  forecast 
sites  based  on  digital  satellite  data.  This  thesis  illus- 
trates the  detailed  analysis  of  the  digital  satellite  infor- 
mation that  can  be  received  by  the  user  in  a  short  (1-2 
hours)  time  frame.  Specifically,  this  work  was  directed  to 
using  the  SPADS  com  cuter  system  toward  process  satellite 
visual  and  infrared  digital  data  and  to  produce  cloud  and 
precipitation   information.      As   SPADS    is      placed    in   the    field 
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forecast  sites,  this  information  will,  for  the  first  time, 
be  made  available  in  an  operational  setting  instead  of  being 
available  only  to  meteorologists  at  the  major  satellite  cen- 
ters. In  the  future,  this  and  other  applications  will  sub- 
stantially  aid  the  operational  forecaster. 
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APPENDIX  k 
SPADS  CLOUD  MODEL  COMPUTER  PROGRAM 


C     CLOD  —  IS  DRIVER/MAIN  PROGRAM 

C      THIS  PROGRAM  ANALYZES  VIS  AND  IR  TOGETHER  TO  DERIVE  TYPES  OP 

C      CLOUDS, CLOUD  HEIGHTS, CLOUD  AMOUNT, CLOUD  TOP  TEMPERATURE, 

C     AND  RELATIVE  INTENSITY  OF  RAIN. 

C     PROGRAMMED  3Y  LT.  C.  NELSON 

C      INITIALIZATION 

C 

COMMON/ A/ ML ( 128), CIRRUS, K,K1 , I SEL 

COMMON/ B/CTYPE( 64) ,PRCIPI(64) ,MSAVG(64) ,M 

COMMON/C/SIGMA(64),SIG1,SIG2,SIG3 

COMMON/ D/ HE IGHT (64), CTEMP ( 64 ) , ACDAVG( 64 ) , E 1 , E2 , ALAT , ALON , J 1 , 
*TLV(11),E3 

COMMON/ F / ZP ,XCEN, YCEN , FLIP , Al , BASLOM , LVL( 1 1 ) 

COMMON/ G/STDP( 11) 

REAL  ML , MS2 ( 8,8), MSAVG , SIGMA , E 1 , E2 , E  , ACLD , ACDAVG , AMCLD , SUM , SUMSIG , 
* HEIGHT, N, CTEMP, SIG1 ,SIG2,SIG3 

INTEGER  CTYPE, CIRRUS, PRCIPI ,M,0,? ,L,PCLD1 ,PCLD2,Q,ISEL,MS(256,3) , 

*VISFN(7) ,IRFN(7) ,MS1(8) 

r 

11=1 

ISEL=0 

M=0 

N=0.0 

El-0.9 

E2=0.55 

E3-O.0 

SUM«0.0 

SUMSIOO.O 

AMCLD*O.0 

CIRRUS-0 

ACLD-0.0 

P«  64 

SIG1=2.5 

3102=20.0 

SIG3-1.5 

PCLD 1=20.0 

0=512 
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c 

C     MENU 
C 

TYPE  "<15>" 

TYPE  "THIS  PROGRAM  IDENTIFIES  CLOUD  TYPES" 

TYPE  "ADDITIONAL  CLOUD  DATA  CAN  ALSO  <JE  OBTAINED" 

TYPE  "THESE  ARE  -  " 

TYPE  "0  NO  ADDITIONAL  INFORMATION1* 

TYPE  "1  PRECIPITATION  INTENSITY" 

TYPE  "2  CLOUD  TOP  TEMPERATURE" 

TYPE  "3  CLOUD  TOP  HEIGHT" 

TYPE  '"4  ALL  OF  ""HE  ABOVE" 

•^vpi7  "q  j?ytt  PROGRAM" 

ACCEPT  "ENTEP  REOUEST31  "  ISEL 

IF(ISEL.EQ.8)GO  TO  30 

IF( ISEL. EQ.O. OR. ISEL. EQ.1)G0  TO  181 

CALL  0PEN(4,"TLVN.0",1,IER) 

IF(IER.E0.1)G0  TO  170 

TYPE  "TLV  FILE  OPEN  ERROR,  RC=  ",IER 
1 70  CALL  RDBLK( 4 , 1 , MS , 1 ,IER) 

IF(IER.EQ.l)GO  TO  ISO 

TYPE   "TLV    FILE    READ    ERROR,    RC=    " , IE" 

CO    TO    70 

no  00  1  1=1,11 

TLV( I )=FLOAT( MS( I ) ) / 100 . 0 

1  CONTINUE 

WRITE<  12 ,826) ( TLvv'(L) , L=  1 , 1 1 ) 
826    F0RMAT(2X,11F6.2) 

2  F0RMAT(S13) 
181    TYPE   "<15>" 

ACCEPT   "ENTER    IR    FILENAME" 
READ(11,2)IRFN(1) 
CALL    0PEN(2,IREN,1 ,IFR) 
IF(IER.EQ.1)G0   TO    190 

TYPE   "IR    FILE   -   OPEN    ERROR,    PC=    " , IER 
ACCEPT    "TRY    AGAIN?", IX 
IF(IX.EQ.1)G0   TO    181 
IF(IX.EO.O)GO   TO    70 
1^0   ACCEPT   "ENTER    VIS    FILENAME" 
READ(11,2)VISFN(1) 
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CALL  0PEN(1,VISFN,1,IER) 

IF(IER.EQ.1)G0  TO  200 

TYPE  "VIS  FILE  -  OPEN  ERROR,  RC=  ",IER 

ACCEPT  "TRY  AGAIN?", IX 

IF(IX.EO.l)C0  TO  190 

GO  TO  69 
C 

C     OPEN  FILES  FOR  STORAGE  OF  CALCULATED  VALUES 
C 

200  CALL  OPEN(5,"CALC1D.O",3,ISR) 
IF(IER.S0.1)CO  TO  201 

TY°E  "DATA  1  FILE  OP^I  ERROR,  RC»  " , IER 
GO   TO  60 

201  CALL  OPEN(7,"CALC2D.0",3,IER) 
IF(IER.S0.1)GO  TO  202 

TYPE  "DATA  2  FILE  OPEN  ERROR,  RC=  " , IER 

GO  TO  69 
C 

C      00  LOOP  TO  CALL  IN  SERIES  OF  DATA 
C     BLOCKS  RUN  FROM.  0  TO  511 

C  LOOP  STARTS  AT  102,192  FOR  IR  TO  MATCH  VIS  (COLLOCATED  CENTER 
C  POINT).  IR  IS  2  MILE  THEREFORE  NEED  TO  USE  EVERY  OTHER  VALUE 
C 

202  13=0 
12=0 

DO  32  11=1,120,2 

I3=(I1-1)+191 

12=1+12 

DO  6  1=1, P 

HEIGHT(I)=0.0 

PRCIPI(I)=0.0 

ACDAVG(I)=O.D 

SIGMA(I)=0.0 

CTYPE(I)=0 

CTE^(I)=0.0 
5  CONTINUE 
C 

C     DRIVER  MODULE  WHICH  WILL  PREPARE  VISUAL  DATA  TO  IE  USED  IN  OTHER 
C     MODULES.   THE  AVERAGE  VISUAL  BRIGHTNESS,  THE  STANDARD  DEVIATION  OF 
C     VISUAL  BRIGHTNESS,  AND  THE  CLOUD  AMONT  IN  BX8  GRID  SETS  'ILL  HE 


C     DETERMINED. 

C      IR  DIGITAL  VALUE  FILE 

C 

CALL  RDBLK(2,I3,MS(1,1) ,1,IER) 

TF(IER.E0.1)C0  TO  203 

TYPE  "IR  FILE  -  READ  ERROR,  RC=  ",IER 

GO  TO  69 

203  GALL  UPK(MS(96,1),MS(1,2),128) 
00  100  1-1,123 

ML( I) -FLOAT (MS (1,2)) 
100  CONTINUE 
Kl-0 

Ji-(i2-i)*a 

GALL  RD5LK( 1 ,J1 ,MS,8 , IER) 
IF(IER.E0.1)G0  TO  204 

TYpE  "visual  FILE  READ  ERROR,  RC=  ",ISR 
GO  TO  69 

204  CONTINUE 

)0   30   £2-1,128,2 

Kl-Kl+1 

K-K2+1 
C 

G  THIS    LOOP    PROCESSES    VISUAL    DATA    fJY    SXS    GRID 

C 

DO   8    1-1,8 

MSKD-0.0 

DO   7   J2-1.8 

V!S2(I,J2)=0.0 
7  CONTINUE 
3  CONTINUE 

L1=(K1-1)*4+1 

DO  11  0=1,3 

CALL  UPK(MS(L1,Q),MS1,3) 

00  10  L»l,8 

MS2(L,0)  =  FL0AT(.!S1(L)) 

10  CONTINUE 

11  CONTINUE 
Nl-8 

do  n  L-1,3 
DO  12  0-1,  ?ll 
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N-N+1.0 

SUM=SUM+MS2(L,Q) 

SUMSIG»SUMSIG+OlS2(L,Q))**2.0 

IF(MS2(L,Q).GE.PCLD1)ACLD«1.0 

IF(MS2(L,0).LT.PCLD1)ACLD=0.0 
AMCLD=ACLD+AMCLD 

12  CONTINUE 

13  CONTINUE 

SIGMA(  Kl  )  =  SQRT(  (  N*SUMSIG-SUM**  2  . 0 )  /  (  M*  (  N- 1 . 0 ) ) ) 

MSAVG(K1)=SUM/N 

CALL  CLOUD 

ACnAVG(Ki)=AICLO/N 

IF<CTYPE(Kl).EQ.8.0R.CTYPE(Kl).Eq. ?) £3=1.0 

IF<ISEL.EO.O.OR.ISEL.EQ.1)GO   TO    18 

CALL   TEMHT 

IF(ISEL.EQ.8)GO   TO    69 
18  CIRRUS-G 

N*0.0 

ACLD-0.0 

AMCLD-0 . 0 

SUM*0.0 

SUMSIG-0.0 

E3=0.0 
30   CONTINUE 

C  STORE    ALL   VALUES 

C  VALUES    ARE    STORED    IN   CALC1    FILE   —   CTYP^RCIPI.CTE 

C  HEIGMT( 1-256   WD) 

C     VALUES  IN  CALC2  —  MSAVG, SIGMA, AC!HVG,0( 1-256  WD) 

C 

DO  23  1=1,64 

MS(I,1)=CTY?E(I) 

J=I+64 

MS{J,1)*PRCIPI(I) 

K5=J+64 

MS(K5,l)»IFIX(CTEMP(I)+.5) 

K4=K5+64 

MS(K4,l)*IFIX(HEIGHT(I)+.5) 

MS(I,2)=IFIX(MSAVO(I)+.5) 
MS(J,2)«IFIX(100.0*SIGMA(I)+.5) 
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23 


22 


26 


MS(K5,2)=IFIX(100.0*ACDAVO(I)+.5) 

MS(K4,2)=0 

CONTINUE 

K7-I2-1 

CALL  WRBLK( 5 ,K7 ,MS( 1 , 1 ) , 1 , IER) 

IF(IER.EQ.l)GO  TO  22 

TYPE  "WRITE  ERROR  5  RE=  ",IER 

GO  TO  69 

CALL  WRBLK(7,K7,MS(1,2) ,1,IER) 

IF(IER.E0.1)GO  TO  26 

TYPE  "WRITE  ERROR  7  RC=  " , IER 

GO  TO  60 

WRITE(12,27)(HEICHT(I),I=l,64) 

WRITE ( 12 ,24)(SIGMA(I) ,1=1 ,64) 

WRITE( 12, 27) (CTEMP(I), 1-1,64) 

WRITE(12,24)(MSAVG(I),I=1,64) 

WRITE(12,24)(ACDAVG(I),I=l,64) 

-/RITE(12,25)(CTY°E(I)  ,1=1 ,64) 

WRITE(12,25)(PRCIPI(I),I=1,64) 

WRITE( 12, 27 )(ML( I), 1-2,123,2) 

24  FORMAT(/,4(/,16(lX,F7.4))) 

25  F0RMAT(/,64I2) 

27  FORMAT(/,4(/,16(lX,F7.2))) 
ACCEPT  "CONTINUE?", ISTP 
IF(ISTP.EQ.1)G0  TO  69 
CONTINUE 

CALL  CLOSE (1 ,IER) 
CALL  CL0SE(2,IER) 
GAEL  CLOSE (4, IER) 
CALL  CLOSE (5, IER) 
CALL  CLOSE (7, IER) 
STOP 
END 


32 
69 
70 


30 


127 


SUBROUTINE  CLOUD 
C 

C     CLOUD  TYPING  MODULE,  ALSO  THE  PRECIPITATION  MODULE  IS  CALLED 
C     FROM  HERE  WHEN  PRECIPITATION  CLOUDS  ARE  IDENTIFIED. 
C 

COMMON/ A/ML ( 128), CIRRUS, K,K1,ISEL 

COMMON/ B/CTYPE( 64) ,PRCIPI(64) ,MSAVG(64) ,M 

REAL  ML , MSAVG ,  HEIGHT ,  CTEMP ,  S IG ,  S IG 1 ,  SIG2 ,  S IG.3 

INTEGER  I , CTY°E ,PRCIPI , CIRRUS , ISEL 
C 

C     THESE  ARE  LILJAS  CONVERTED  GOES  THRESHOLDS 
C 

H2=  84.0 

H3-100.0 

H4-118.0 

115  =  15°). 0  • 

H6=255.0 
C 

C     THE  FOLLOWING  IF  STATEMENTS  CHECK  BOTH  IR  (ML)  AMD  VIS  (MSAVG) 
C     VALUES  AGAINST  THRESHOLDS  TO  DETERMINE  CLOUD  TYPE.   THE 
C      STANDARD  DEVIATION  (SIG)  OF  THE  8X8  VISUAL  GRID  IS  USED. 
C 
C     TYPE  1  AND  2  CIRRUS /ALTOSTRATUS 

IF(ML<K).GT.H4.AND.ML(K).LE.H6.AND.MSAVG(K1).GT.5.0.AND. 
*MSAVG(K1).LE.23.0)GO  TO  70O 
C     TYPE  3 FOG/ STRATUS /CUMULUS  HUJIULUS 

IF(ML(K).GT.H2.AND.ML(K).LE.H4.AND.MSAVG(K1).GT. 13.0. AND. 
*MSAVG(K1).LE.20.0)G0  TO  Q00 
C     TYPE  4  AND  5  THICK  FOG/STRATOCUMULUS/CUMULUS  fffJMULUS 

IF(ML( K) .GT.H3 . AND.ML(K) .LE. H4 . AND ,MSAVG( Kl ) .GT. 20 .0 .AND. 
*MSAVG(K1).LE.36.0)GO  TO  600 

GO  TO  601 

600  M=l 

GO  TO  800 
C     TYPE  6  AND  7  SMALL/LARGE  CUMULUS  CONGESTUS 

601  IF(ML(K).GT.H4. AND. ML(K).LE.H5. AND. MSAVG(K1).GT. 23.0. AND. 
*MSAVG(K1).LE.39.0)GO  TO  605 

GO  TO  606 
605  M=2 

GO  TO  800 


C     TYPE  8 NIMBOSTRATUS  (MULTI-LAYERED) 

606  IF(ML(K).GT.H5. AND. ML(K).LE.H6.AND.MSAVG(K1).GT. 23.0. AND, 

*MSAVC(K1).LE.31.0)GO  TO  910 

C     TYPE  9 CUMULONIMBUS 

IF(ML(K). GT.H5. AND. ML(K).LE.H6. AND. JISAVG(K1).GT. 11.0) 
*CO  TO  <U5 
CTYPE(K1)=0 
GO  TO  38 
C     CALL  STDDEV  MODULE  TO  FURTHER  DISCRIMINATE  CLOUD  T'^ES 
700  M=3 

800  CALL  STDDEV 
GO  TO  38 
C     LABEL  TYPES 
300  M=4 

CO  TO  800 
38  PRCIPI(K1)=0 
GO  TO  40 
910  CTYPE(K1)=S 

GO  TO  920 
915  CTYPE(K1)=9 
C      30TN  8  AND  ^  TYPES  ARE  CONSIDERED  RAIN  CLOUDS,  SO  CALL 
C     PRSCIP  MODULE. 

920  IF(ISEL.EQ.0.0R.ISEL.EQ.2.OR.ISEL.EQ.3)G0  ~o  40 
CALL  PRECIP 
40  RETURN 
END 
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SUBROUTINE  STDDEV 
C 

C     STANDARD  DEVIATION  TEST  MODULE  HELPS  DISCRIMINATE  FURTHER  THE 
C     TYPE  OF  CLOUDS.   ALSO  IT  SETS  CIRRUS  FLAG  FOR  USE  IN  THE  TEMHT 
C     MODULE . 
C 

COMMON/ A/ML( 128), CIRRUS, K.Ki.ISEL 

COMMON/B/CTYPE(fi4) ,PRCIPI(64) ,MSAVO( 64) ,M 

COMMON/C/SIGMA<64),SIGl,SIG2,SIG3 

REAL  ' 1L , MSAVG , HEIGHT , CTEMP , SIG , Y, SIG1 , SIG2 , S IG3 

INTEGER  CTYPE,PRCIPI , CIRRUS , I ,M 
C 

IF(M.EQ.4)G0  TO  818 

IF(M.EQ.l)GO  TO  8i0 

IF(M.E0.2)CO  TO  313 
C     CIRRUS  VS  ALTOSTRATUS,  FIRST  TEST 
C     THIS  EQUATION  IS  FOR  LILJAS  CONVERTED  GOES  THRESHOLDS 

Y-158. 094-6. 6076*MSAVG(K1)+4.159*MRAVG(K1)**2.0 

IF(Y.LF.ML(K))GO  TO  814 
C      SECOND  TEST  FOR  CI/AS 
C     THIS  EQUATION  IS  FOR  LILJAS  CONVERTED  GOES  THRESHOLDS 

Y»115.24+0.5246*MSAVG(K1)+0.0395*MSAVG(K1)**2.0 

IF(Y.GE.ML(K))GO   TO   816 
C  THIRD   TEST    FOR   CI/AS 

IF(SIGM.A(K1);GT.SIG3)G0   TO   814 

CO   TO   816 

810  IF(SIGMA(K1).LT.SIG1)G0  TO  811 
C     TYPE  5  —  CUMULUS  HUMULUS 

CTYPE(Kl)-5 
GO  TO  820 
C     TY^E  4  —  THICK  FOG/ STRATO CUMULUS 

811  CTYPE(K1)=4 
GO  TO  320 

313  IF(SIOMA(Kl).LT.3IG2)GO  TO  817 
CTYPE(K1)«6 

C     TYPE  6  —  SMALL  CUMULUS  CONGESTUS 

GO  TO  820 
C      TYPE  1  —  CIRRUS 

314  CTY?E(K1)=1 
CIRRUS=1 


GO  TO  820 

TYPE  2  —  ALTO STRATUS 

816  CTYPE(K1)=2 
00  TO  820 

TYPE  7  —  CUMULUS  COMGESTUS  -  LARGE 

817  CTYPE(K1)=7 
GO  TO  820 

818  IF(SIGMA(K1).LT.SIG1)G0  TO  819 
TYUE  5  —  CUMULUS  HUMULUS 
CTYPE(K1)=5 

GO  TO  820 

TYPE  3  —  FOG/ STRATUS 

819  CTYPE(K1)«3 

820  RETURN 
END 
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SUBROUTINE  PRECIP 
C 

C     THIS  MODULE  IDENTIFIES  PRECIPITATION  INTENSITY  AREAS 

C      1-LIGHT,  2=MODERATE,  3=HEAVY,  AND  0«NO  RAIN. 

C 

COMMON/ A/ML( 128), CIRRUS, K,K1,ISEL 

COMMON/ 3/CTYPE( 64 ) ,PRCTPI( 64)  ,MSAVG(64)  ,M 

REAL  ML ,  HE  I GHT ,  CTE  MP ,  MS AVG 

INTEGER  K, CIRRUS , PRCIPT , CTYPE , SUM1 
C 

C     THESE  ARE  LILJAS  CONVERTED  GOES  THRESHOLDS 
C 

11=134 

12=195 

13-224 
C 

C     ADD  IR  (ML)  AND  AVERAGE  VISUAL  (MSAVG) ,  THIS  SUM  WILL  BE  TESTED 
C      FOR  INTENSITY  OF  RAIN. 

SUM1»ML(K)+MSAVG(K1) 

IF(SUM1.GE. II. AND. SUM1.LE. 12)00  TO  43 
IF(SUM1.0T.I2.AND.3UM1.LE.I3)G0  TO  44 
IF(SUM1.GT.I3)G0  TO  45 
PRCIPI(K1)=0 
GO  TO  49 

43  PRCIPI(K1)=1 
GO  TO  49 

44  PRCIPI(K1)»2 
GO  TO  49 

45  PRCIPI(K1)=3 
49  RETURN 

END 
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SUBROUTINE  TEMHT 
C 

C  THIS  SUBROUTINE  CALCULATES  THE  CLOUD  TO?  TEMPERATURE  (CTEMP) 

C  USING  THE  IR  RADIANCE  FOR  THE  SURFACE  (NCLR),  THE  AMOUNT  OF  CLOUD 

C  (ACDAVG),  EMISSIVITY  (E),  AND  THE  MEASURED  IR  RADIANCE  FROM  THE 

C  SATELLITE  (ML).   CTEMP  IS  THEN  CORRELATED  TO  AN  AREA 

C  REPRESENTATIVE  UPPER  AIR  SOUNDING  (UALV)  TO  DETERMINE  THE  CLOUD 

C  TOP  HEIGHT  (HEIGHT). 

COMMON/ A/ ML( 128), CIRRUS, K,K1,ISEL 

COMMON/ D/HEICHT(  64)  ,CTEMP(64)  ,ACI)AVG(64)  ,E1  ,E2  ,ALAT,ALON,Jl  , 
*TLV(11) ,E3 
COMMON/ G/STDP( 11) 

REAL  : !L , HEIGHT , CTEMP , £ 1 , E2 , NCLR,NCLD , ACDAVG , E , SIG 1 , SIG2 , SIG3 , 
*TLV,STDP 
INTEGER  I, CTYPE,PRCI?I, CIRRUS 

DATA  STDP/ 1010.0, 1000 .0,850.0,700.0, 500 . 0 , 400 . 0 , 300 . 0  ,250 . 0 , 200 . 0  , 
*150. 0,100.0/ 
E=0.0 
C     CHECK  TO  SEE  IF  THERE  ARE  ANY  CLOUDS,  IF  NOT  SKIP  TO  END 

IF(ACDAVG(K1).LT.0.1)GO  TO  61 
C     CONVERT  SURFACE  TEMPERATURE  (TSFC)  TO  NCLR 
IF(TLV(1).CE.242.0)NCLR=2.0*(331.0-TLV(1)) 
IF(TLV(1).LT.242.0)NCLR=420.0-TLV(1) 
C      IS  CIRRUS  FLAG  SET?   IF  YES,  USE  E2  FOR  E.   IF  NOT,  USE  El  FOR  E. 

IF( CIRRUS. EO.l)G0  TO  850 
C      IS  THE  CB  OR  NS  CLOUD  FLAG  SET?  IF  YES  USE  E3  FOR  El. 
IF( S3. EO. 1.0)GO  TO  849 
E-El 

GO  TO  851 
849  E=E3 

GO  TO  851 
350  E=E2 
C     COMPUTE  CLOUD  'RADIANCE 

851  NCLD»ML(K)-NCLR/(E*ACDAVG(K1))+NCLR 
C      CONVERT  TO  CTEMP 

IF( NCLD.LE. 1 78 .0)CTEMP( Kl )=33 1 . 0-NCLD*0 . 5 
IF(NCLD.GT.173.0)CTEMP('<1)=420.0-NCLD 
WRITE(12,353)NCLR,NCLD 
353  F0RMAT(/1X,5(1H*,E15.3),"NCLR,NCLD") 
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IF(ISEL.EQ.2)G0  TO  60 
C 

C     INTERPOLATE  UPPER  AIR  SOUNDIMG  TO  CONTINUOUS  CURVE  OF  TEMP 
C     VERSUS  HEIGHT,  THEN  COMPARE  TO  CTEMP  TILL  GET  MATCH.   CHECK 
C     FROM  TOP  DOWN.  -THEN  USE  MATCH  TO  SELECT  CORRECT  HEIGHT  FOR 
C     CLOUD  TOP. 
C 

DO  312  13=2,11 

IA-13-I3 

IF(IA.E0.11.AND.TLV(11).OT.CTEMP(K1))00  TO  316 

IF(TLV(IA).LT.CTEMP(fU))CO  TO  313 

T~(TLV(IA).EO.CTEMD('U))C-0  TO  314 

312  CONTINUE 

313  IF(IA.E0.2)IA=3 

RHG=EXP(  ALOG(STDP(  IA+1 )  )-(  ALOC(  ST">?(  IA+1 )  )-ALOC(  ST^P(  TA)  )  )  * 
*(TLV( IA+1) -CTEMP (Kl ) ) / ( TLV( IA+1 )-TLV( IA) ) ) 
"GO  TO  315 
816  TYPE  "TEMPERATURE  IS  ABOVE  100  MB" 
ACCEPT  "CONTINUE?", ICOMT 
IF(ICOMT.EO.O)ISEL=8 

314  RHC=STDP(IA) 

815  HEIGHT(Kl)»RHG+0.5 

60  RETURN 

61  CTEMP (K1)=TLV(1) 

HEIGHT(K1)»STDP(1)' 

RETURN 

END 
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